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SUMMARY 

Several computer subroutines have been developed that Interpolate three 
types of nonanalytlc functions: univariate, bivariate, and map. The routines 

use data In floating-point form. However, because they are written for use on 
a 1 6-b 1 t Intel 8086 system with an 8087 mathematical coprocessor, they execute 
as fast as routines using data In scaled Integer form. Although all of the 
routines are written In assembly language, they have been Implemented In a 
modular fashion so as to facilitate their use with high-level languages. 


INTRODUCTION 

The use of digital microcomputers In real-time simulation and control 
applications has created a need for high-speed function generation routines. 
These routines are needed to mathematically simulate the relationships between 
the Inputs and outputs of a physical system. 

For function generation routines being used In the control of systems 
such as gas turbine engines, speed of execution Is a major concern. Previ- 
ously, because of this speed requirement, function generation programs could 
only deal with data In an Integer format (ref. 1). However, by using a mathe- 
matical floating-point hardware coprocessor, these routines can be extended to 
deal with data In floating-point form. 

The programs described In this report are written In assembly language 
for use with the Intel 8086 microprocessor and the 8087 coprocessor chip. 

They use straight-line Interpolation as the numerical estimation method and 
are callable from high-level languages. The ability of these programs to be 
directly linked with a high-level language makes It possible to develop the 
main logic of a control In a high-level language with calls to the assembly 
language routines when high-speed function generation Is needed. 

These programs deal with three types of functional relationships: (1) 

simple Input-output functions (univariate functions), (2) functions of two 
variables with Identical x data values on each curve (bivariate functions), 
and (3) functions of two variables with different x data values on each curve 
(map functions) . Figures 1 to 3 Illustrate each type. The programs are 
described by examining the three types of functions, by discussing the computer 
subroutines written to Implement these functions, and by examining how the 
high computation speeds were achieved. 





DESCRIPTION OF FUNCTION TYPES 


The six programs described In this report can Interpolate outputs for 
three kinds of functional relationships. These functional relationships are 
described In this section. 


Function Type One (Univariate Function) 

This Is the simplest type of function. It Is a simple Input-output 
relationship In which all data points are on the same curve (fig. 1). Every 
Input value x uniquely Identifies an output value y. The function can be 
symbolized as Y y = F(Xy). An output value Y y for which there Is no corre- 
sponding experimental data point Xy can be found by linearly Interpolating 
between adjacent data points. This Interpolation Is carried out by using 
equation (1). 



The relationships between variables are shown In figure 1. 


Function Type Two (Bivariate Function) 

This function consists of a family of curves (fig. 2) Instead of the 
single curve of function type one. Each curve has a particular value of z 
assigned to It. Therefore It Is necessary to provide two Input values, Xy and 
Zy, before an output can be determined. This type of functional relationship 
can be symbolized as Yy = F(Xy, Zy) . For this fun ctio n a restriction Is 
made that each curve must have the same number of data points. Additionally, 

It Is required that the data points for each curve be taken at Identical values 
of x. The output Yy Is found by linearly Interpolating between adjacent 
x and z values by using equations (2) to (4): 



The relationships between variables are shown In figure 4. 
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Function Type Three (Map Function) 

This last function, which also consists of a family of curves, Is the most 
general one. The restriction requiring all curves to have the same x break- 
points Is removed. However, each curve Is still required to have the same 
number of breakpoints per z curve (fig. 3). Because each curve has a unique 
set of x and y breakpoints, additional calculations must be made. The 
equations required to obtain the Interpolated output Yy = F(Xy, Zy) are 



The relationships between variables are shown In figure 5. 


DESCRIPTION OF COMPUTER ROUTINES 

The programs used to Implement function generation were developed for the 
Intel 8086 1 6— bl t microprocessor coupled with an 8087 mathematical coprocessor. 
The 8086 chip uses 16-bit address registers and 16-bit segment registers. An 
address register and a segment register are combined to form a 20-bit address 
that can access up to 1 megabyte of storage (fig. 6). In Intel "small-memory 
model" programs the segment registers remain constant while the program Is 
executing. Although this restriction allows the program to address only 64K 
bytes of storage, It simplifies many programming tasks and has the least amount 
of programming overhead. Intel "large-memory model" programs do not require 
the segment registers to remain constant. Thus large-model programs must set 
up and maintain the segment registers at their correct values. Because of the 
subtle differences between the large and small models, six separate function 
generation routines are presented In this report, a large- and small-model 
program for each of the three types of function. Further Information on the 
large- and small -memory models for the 8086/8087 processors Is given In 
references 2 to 5. 
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Linking to High-Level Languages 

The function generation routines can be linked to both large- and small- 
model high-level language programs because all Intel conventions regarding the 
saving and restoring of registers are followed. This means that the routines 
can be used with Fortran-86, a large-model, high-level language, or PL/M-86, 
which can be either a large- or small-model language. 

The high-level language must supply the function generation routines with 
two types of Input parameter: (1) the address In memory of a data block hold- 

ing Information about the experimental breakpoints (the breakpoint Information 
block), and (2) the Input values of the function. 

The breakpoint Information block contains two types of Information: (1) 

Information used to describe the format of the experiment data, and (2) the 
arrays holding the data. All Items In the breakpoint Information block are 
Identified In figure 7. 

An example of a typical high-level-language call Including setup of the 
breakpoint Information block Is shown In figure 8. According to parameter 
passing conventions, the first storage address In the breakpoint Information 
block (the variable XPTR) Is passed by the high-level-language program on the 
8086 stack, and the real-number Input values are passed on the 8087 stack. 

The final Interpolated result produced by the lookup routine Is returned 
to the high-level-language program on the 8087 stack. When the parameters are 
passed to the subroutine as shown In figure 8, this process Is taken care of 
automatically by the compiler and assembly language function routines and Is 
transparent to the high-level language programmer. 


Out-Of-Range Conditions 

For a typical control application an out-of-range data Input must not be 
allowed to cause an unpredictable output. This Is necessary to prevent damage 
to the control system that could occur If, for Instance, an actuator were 
pushed beyond Its range of operation. If a value Is out or range, either too 
high or too low, the programs will output the closest boundary value of the 
function. If, for example, the Input value Is higher than the last experi- 
mental data point In the breakpoint Information block of a simple Input-output 
curve, the routine will output the highest experimental value. If the Input 
z value Into a two-variable function Interpolation routine Is above the 
maximum z data value, the routine will use the highest allowable z curve 
when performing the Interpolation (fig. 9). 

Dealing with out-of-range values In this way keeps the output from 
becoming unpredictable and results In a control system that Is well behaved. 


OPERATION OF COMPUTER ROUTINES 


Six subroutines are described In this report, two for each type of 
function: They are RFUN1S, RFUN1L, RFUN2S, RFUN2L, RFUN3S, and RFUN3L. The 

last letter In the program name tells the type of memory model with which the 
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program will Interface. A "S" suffix Indicates that the routine Is to be used 
In the small-memory model, and a "L" Indicates use with the large-memory model. 
Listings of all programs, Including example calls, are Included In the appendix 
to this report. 


RFUN1S and RFUN1L Programs 

The routines RFUN1S and RFUN1L, which Interfaces with the Intel small- and 
large-model programs respectively, deal with simple univariate Input-output 
function lookups. The data are set up In the breakpoint Information block as 

shown In figure 10. Figure 8 shows how these block data are set up by a 

Fortran calling program. 

XPTR Is an Index pointing to the location of the data In the most recent 
previous call to the routine. ZPTR Is not used In these routines but Is 

Included to make the data structure similar to that used In RFUN2 and RFUN3. 

NXPTS Is the number of data breakpoints, NZPTS Is not used but Is Included for 
the same reason as ZPTR. XARRAY, which contains the breakpoint values of x, 

Is stored next followed by YARRAY, the function output for each x value. 

During execution three basic operations are performed. First the routine 
Initializes Its registers. Then It determines whether the Input value Is 
higher or lower than the Initial x value Indicated by the XPTR Index. The 
program then scans the XARRAY table to find the x segment that corresponds to 
the x Input value. Once the scanning Is complete, the 8087 chip Is used to 
perform the Interpolation by using equation (1), which was discussed earlier. 
Finally, the routine updates the Index pointer and returns to the calling 
program with the Interpolated function output value y on the 8087 stack. 


RFUN2S and RFUN2L Programs 

RFUN2S and RFUN2L generate values for two-variable functions with Identi- 
cal breakpoints (as In fig. 2) In the small- and large-model programs, respec- 
tively. They function similarly to RFUN1S and RFUN1L with one additional step. 
The routines must perform a search through the z values as well as through 
the x values. They must also read and update the ZPTR value In the break- 
point Information block (this variable was not used In the RFUN1 programs). 

The data are stored In this program as shown In figure 11. Since all of 
the x breakpoints occur at Identical points on all curves, only a one- 
dimensional XARRAY Is necessary. However, because each curve will have a dif- 
ferent y output value at each breakpoint, a two-dimensional YARRAY Is needed. 
The arrangement of this array Is such that the first row corresponds to the 
outputs of the first z curve, the second row corresponds to the outputs of 
the second z curve, and so on, up to the number of curves that exist. 

The programs make use of equations (2) to (4) to compute the proper Inter- 
polated output. As In the RFUN1 programs, the final Interpolated function 
output Is returned on the 8087 stack. 
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RFUN3S and RFUN3L Programs 


The RFUN3 programs generate outputs for two-variable functions with Indi- 
vidual breakpoints (as In fig. 3). They operate similarly to RFUN2S and RFUN2L 
except that additional calculations are made during each x scan to find the 
breakpoints on a new z curve. The RFUN3 programs first scan the ZARRAY to 
determine which curves the z Input value Is between. Next It uses 
equations (5) and (6) to generate the x breakpoints of this new Interpolated 
z curve. The x value search proceeds by generating new x values as a 
function of z until the upper and lower boundaries of the Input x are 
found. Then Interpolation between the generated segment endpoints by using 
equations (7) to (9) determines the final output value y. Finally the x and 
z Index pointers are updated and the y output value Is returned to the 
calling program on the 8087 stack. 

The data setup for the RFUN3 programs Is more complicated than that for 
either the RFUN1 programs or the RFUN2 programs because each breakpoint Is 
Independent of all others. Not only does each curve have an Individual set of 
y output values, as In the RFUN2 programs, but each also has an Individual set 
of x values. To support this arrangement, the data are set up as shown In 
figure 12. XARRAY Is a two-dimensional array containing the breakpoints for 
each curve. The first row corresponds to the first z curve, the second row 
to the second z curve, and so on as In the YARRAY used In the RFUN2 programs. 
The YARRAY used In the RFUN3 programs Is set up the same as the YARRAY In the 
RFUN2 programs. 


IMPLEMENTING FAST EXECUTION SPEEDS 

Since function generation routines are used most often In real-time appli- 
cations where execution time Is critical, an efficient algorithm Is very 
Important. Several techniques are used In the function generation routines to 
attain this speed efficiency. 

Before Interpolation Is possible, the function generation routine must 
determine the magnitude of the Input by searching through the breakpoint 
Information block values for the one closest to the Input value. Typically, 
since a specific function routine Is called only once per control update 
Interval, and since the Input value normally changes relatively slowly from one 
control Interval to the next, the breakpoint Information block search Is sped 
up If the Index location of the previous search Is known. All of the function 
generation programs described In this report store the final Index value for 
use as the starting point of the next search. This Index value, XPTR for 
example, Is stored at the beginning of the breakpoint Information block and Is 
updated at the end of each call of the function generation routine. 

Another way to decrease the execution time of the routine Is to use memory 
access Instructions only when necessary since these take longer to execute than 
register access Instructions. The routines described In this report store any 
frequently used Integer value .In an 8086 Internal register and any frequently 
used real value In an 8087 register, thus eliminating unnecessary memory 
accesses. 
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A third approach used to decrease the execution time Is to use In-line 
coding of Instructions rather than subroutine calls. In-line coding may cause 
the size of the code to be larger, but It will execute faster because the 
processor does not have to store pointers and status registers as It does for 
subroutine calls. 


Function Routine Execution Times 

Typical execution times for the entire lookup routine were experimentally 
determined to be 0.50 ms for a univariate function lookup, 0.90 ms for a 
bivariate function lookup, and 1.00 ms for a map function lookup. 


CONCLUDING REMARKS 

Three types of nonanalytlc function routines have been developed for use 
with control algorithms written In high-level languages. These routines 
operate on data In floating-point form and are fast enough to be used In real- 
time control systems. Additionally, problems associated with the 8086 
segmented architecture have been resolved so that the routines can be easily 
and efficiently Integrated Into most control schemes. 
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8086/87/88/186 MACRO ASSEMBLER RFUN1S 


08/30/84 PAGE 1 


SERIES-III 8036/87/88/186 MACRO ASSEMBLER V2.0 ASSEMBLY OF MODULE RFUN1S 
OBJECT MODULE PLACED IN :F1:RFUN1S.0BJ 
ASSEMBLER INVOKED BY: ASN36.86 :FlsRFUNlS.SRC 


LOC OBJ LINE SOURCE 


0000 

0002 

0004 

0006 


1 

2 

3 

4 

5 

6 
7 
3 

9 

10 
11 
12 

13 

14 

15 

16 
17 
IS 

19 

20 
21 
22 

23 

24 

25 

26 
27 
23 

29 

30 

31 

32 

33 

34 

35 

36 

37 
33 

39 

40 

41 

42 

43 

44 

45 

46 

47 
43 

49 

50 




* ROUTINE NAME: RFUN1S * 

* VERSION: 2.0 * 

* DATE: AUG. 9, 1933 * 

* PROGRAMMER'S NAME: MICHAEL MACK IN * 

* PURPOSE: TO INTERPOLATE A SIMPLE * 

* FUNCTION FOR SMALL-MODEL * 

* PROGRAMS * 




THIS SUBROUTINE OUTPUTS THE ESTIMATED VALLE OF A FUNCTION YV = F(XV) 
BY INTERPOLATING FROM KNOWN SOLUTION POINTS? USING THE RELATION 
YV = ((XV-XL)/(XH_XL) ) (YL-YL) + YL WHERE 
YV = ANSWER TO BE COMPUTED 
XV = INPUT VALLE 
XH = KNOWN HIGH X POINT 
XL = KNOWN LOW X POINT 
YH = KNOWN HIGH Y POINT 
YL = KNOWN LOW Y POINT 

INPUT REQUIREMENTS: 

1) THE FOLLOWING ADDRESSES WILL BE PASSED TO THIS SUBROUTINE: 

A. ON TOP OF THE 3036 STACK: 

THE ADDRESS (OFFSET) OF A DATA AREA WITH THE FOLLOWING- 
FORMAT 

DSTRUC STRUG 

XPTR DU 
2PTR DW 
NXPTS DU 
NZPTS DW 
X ARRAY DD 
YARRAY DD 

DSTRUC ENDS 

B. ON TOP OF THE 3037 STACK: 

THE INPUT VALUE? XV? TO BE INTERPOLATED 

2) IT IS REQUIRED THAT THE 8087 CHIP HAVE 3 EMPTY REGISTERS 

WHEN THIS INTERPOLATION ROUTINE IS CALLED. 

OUTPUT EFFECTS: 

1) REGISTERS DESTROYED: AX? BX? CX, BX, DI» SI 

2) INTERPOLATED RESULT RETURNED ON TOP OF 3087 STACK 

3) XPTR LOCATION UPDATED TO INDEX X VALUE JUST PRECEDING XV. 

(TO SPEED FUTURE SEARCHES FOR X) 


? 5X ARRAY INDEX 

? ;Y ARRAY INDEX 

? 5W0. OF ELEMENTS IN X ARRAY 

? ?NG. OF ELEMENTS IN Y ARRAY 

NXPTS DUP(?) ’X ARRAY 

NYPTS DUPf?) !Y ARRAY 
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8086/87/88/186 MACRO ASSEMBLER 


RFUN1S 


08/30/84 PAGE 2 


LOC OBJ LIME SOURCE 


oooo 

0002 

0004 


0000 ???? 


51 

52 

53 

54 

55 

56 

57 

58 
5? 
60 
61 
62 

63 

64 

65 

66 
67 
63 
6? 

70 

71 

72 

73 

74 

75 

76 

77 
73 
7? 
80 
81 
82 
33 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 
93 

99 

100 
101 
102 

103 

104 

105 


EXAMPLE PL/M PROGRAM USING RFUN1S: 

♦SMALL 
PU1PGM: DO; 

DECLARE DSTRUC STRUCTURE ( 

XINDEX INTEGER, ZINDEX INTEGER, NXPTS INTEGER, NZPTS INTEGER, 
XARRAYI7) REAL, YARRAYI7) REAL) PUBLIC 
INITIAL (0,0, 7, 7, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 0.0, 1.0, 
4.0,9.0,16.0,25.0,36.0); 

RFUN1S: PROCEDURE ( X , STRUC.ADDR ) REAL EXTERNAL; 

DECLARE X REAL; 

DECLARE STRUC-ADBR POINTER; 

END RFUN1S; 

PLMTEST 2 PROCEDURE PUBLIC; 

DECLARE I INTEGER; 

DECLARE (X»Y> REAL; 
x = 1.5; 

Y = RFUN1SCX, 6DSTRUO? 

CALL PRINT-ANSWER ( Y > " 

END; 

END FUNPGM; 

END PLMPGM? 

m««mm*M*#*m***™mm****mm************************** 


NAME RFUN1S 
CGROUP GROUP CODE 

CGROUP GROUP DATA 

ASSUME CS: CGROUP, DSsDGROUP 
PUBLIC RFUN1S 

;m**M*#M*mmM**«**M****M*m*#M***Mm*m™*#*****m**# 


STACK STRUCTURE 


PARAMS 

STRUC 



OLD-BP 

DU 

7 

! SAVED BP REGISTER 

RETURN 

DU 

7 

! RETURN ADDRESS 

DADDR 

DU 

7 

! ADDRESS OF DATA STRUCTURE 

PARAMS 

ENDS 



;*m*«**M****m«*M**m***»**m*m*»****m****»#*M)HHHHtt**** 

DATA 

SEGMENT PUBLIC 

'DATA' 


STMSAVE DM ? 
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8036/87/88/186 MACRO ASSEMBLER RFUN1S 
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LOC OBJ 


LINE 

SOURCE 





106 




— 


107 

BATA ENDS 





103 






10? 




110 




111 




— 


112 

CODE SEGMENT PUBLIC 'CODE' 




113 




0000 


114 

NEARCQDE 

PROC NEAR 




115 






116 






117 




0000 


113 

RFUN1S: 





11? 



?SAVE REGISTERS 

0000 55 


120 

PUSH 

BP 


0*001 SBEC 


121 

MOO 

BP, SP 


0003 SB5E04 


122 

MOV 

BX, tBPLBABDR 

? GET START ADDRESS OF DATA STRUCTURE 

0006 8B3F 


123 

MOV 

DI» CBXLXPTR 

;DI = ELEMENT INDEX 

0003 DiE7 


124 

SHL 

DI, 1 


000ft D1E7 


125 

SHL 

01, 1 

;BI = BYTE INDEX = ELEMENT INDEX * 4 

OOOC 3B7704 


126 

MOV 

SI, EBX3.NXPT3 

’SI = HXPTS 

OOOF B1E6 


127 

SHL 

SI. 1 


0011 D1E6 


123 

SHL 

SI, 1 

iSI = NO. OF BYTES=NXPTS * 4 

0013 83C30S 


12? 

ADD 

BX, 8 

;BX = START OF X ARRAY 

0016 3BB3 


130 

MOV 

DX, BX 


0018 03D6 


131 

ADD 

DX, SI 

;DX = START OF Y ARRAY 

001 A 9BD811 


132 +2 

FCOM 

DWORD PTR EBX + DI] 

COMPARE ARRAY ELEMENT WITH XV 

00 ID 9BDD3E0000 

R 

133 +1 

FSTSH 

STHSAVE 


0022 9B 


134 +1 

FWAIT 



0023 A 10000 

R 

135 +1 

MOV 

AX, STHSAVE 


0026 250041 


136 +1 

AND 

AX, 4100H 

MASK-IN COMPARISON BITS 

<>D2? 3D0001 


137 +2 

CMP 

AX, 0100H 


002C 7432 


133 +2 

139 +2 

140 +1 

JE 

LEFOF_$TART 

! EXIT IF XV < ARRAY ELEMENT 



141 



IF X > DATA (ORIGINAL INDEX) 

002E 


142 

RIGHT.OF-START: 


*** THEN BEGIN ##* 

002E 83C704 


143 

ADD 

DI, 4 

INDEX = INDEX + 4 

0031 3BFE 


144 

RSI: CMP 

DI, SI 

WHILE (INDEX <= NO. OF BYTES) 

0033 741E 


145 

JE 

LI ! 




146 



*** DO BEGIN *** 

0035 9BD311 


147 +2 

FCOM 

DWORD PTR [BX + DI] ; 

COMPARE ARRAY ELEMENT WITH XV 

0033 9BBD3E00G0 

R 

143 +1 

FSTSH 

STHSAVE 


0030 9B 


149 +1 

FWAIT 



003 E A 10000 

R 

150 +1 

MOV 

AX, STHSAVE 


0041 250041 


151 +1 

AND 

AX, 4100H ! 

: MASK-IN COMPARISON BITS 

0044 300001 


152 +2 

CMP 

AX, 0100H 


0047 740A 


153 +2 

JE 

LI s 

: EXIT IF XV < ARRAY ELEMENT 

004? 3D0040 


154 +2 

CMP 

AX, 4000H 


0040 7405 


155 +2 

156 +1 

JE 

Li ; 

! EXIT IF XV = ARRAY ELEMENT 

004E 83C704 


157 

ADD 

DI, 4 s 

' INDEX = INDEX + 4 

0051 EBDE 


158 

JMP 

RSI ; 

! *** END (NHILE) *** 

0053 83EE04 


15? 

Li: SUB 

SI, 4 j 


0056 83EF04 


160 

SUB 

DI, 4 ; 
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LOC OBJ 


LIKE 

SOURCE 



005? 3BFE 


161 

CMP 

DIi SI 

IF INDEX > NO. OF BYTES 

0053 7453 


162 

JE 

QUUF-RANGE 

' THEN OUT_OF_ RANGE 

0050 EB3290 


163 

JMP 

FOUND-LCW-INDEX 

ELSE FOUND_LOW- INDEX 



164 



*** END (IF) *** 



165 



IF X < DATA (ORIGINAL INDEX) 

0060 


166 

LEFTJ3F-START: 


**# THEN BEGIN *** 

0060 83EF04 


167 

SUB 

DI» 4 

INDEX = INDEX - 4 

0063 33FF00 


163 

LSI: CMP 

DIi 0 

WHILE (INDEX > 0) 

0066 7C1E 


16? 

JL 

Ml 




170 



*#* DO BEGIN *** 

0068 9BD811 


171 +2 

FCOM 

DWORD PTR CBX + DII 

COMPARE ARRAY ELEMENT WITH XV 

006B 9BDD3E0000 

R 

172 +1 

FSTSW 

STWSAVE 


<X'70 9B 


173 +1 

FWAIT 



0071 A 10000 

R 

174 +1 

MOV 

AX, STWSAVE 


0074 250041 


175 +1 

AND 

AX, 4100H ■ 

; MASK-IN COMPARISON BITS 

0077 3B0000 


176 +2 

CMP 

AX, OOOOH * 

! EXIT IF XV > ARRAY ELEMENT 

007 A 740A 


177 +2 

JE 

MI 


007C 3D0040 


173 +2 

CMP 

AX, 4000H 


007F 7405 


17? +2 

JE 

Ml 

i EXIT IF XV = ARRAY ELEMENT 



130 +1 




0081 33EF04 


131 

SUB 

DI> 4 

! INDEX = INDEX - 4 

0034 EBBD 


132 

JMP 

LSI 

; *** END (WHILE) *** 

0036 83FF00 


183 

His CMP 

DI. 0 

? IF INDEX O INDEX-MIN 

008? 7006 


134 

JGE 

FOUND-LOWLINDEX 

5 THEN FOUND-LOW-INDEX 

008B 83C704 


185 

ADD 

DI, 4 

" 

008E EB2590 


136 

JMP 

OUT-OF-RANGE 

; ELSE OUT-OF-RANGE? 



187 



; w END (IF) *** 



183 




0091 


18? 

FOUND-LOW. INDEX 

: 


0091 


190 

INTERPOLATE: 


?XV ALREADY ON TCP BY PLM PGM, 

0091 9BD901 


191 

FLO 

DWORD PTR CBX+DI3 

; STORE XL ON STACK TOP 

0094 9BBCE9 


192 

FSUB 

ST ( 1 ) , ST 

? (XV— XL) ON STACK BOTTOM 

0097 9BDS6904 


193 

FSUBR 

DWORD PTR [BX+DI+41 

? (XH-XL) ON STACK TOP 

009B 9BDEF9 


194 

FDIV 


? < XV-XL ) / ( XH-XL ) LEFT ON STACK 

009E 37BA 


195 

XCHG 

BX, DX 


OOAO 9BD901 


196 

FLD 

DWORD PTR EBX+DI3 

; STORE YL ON STACK TOP 

00A3 9BD94104 


197 

FLO 

DWORD PTR [BX+OI+4] 

; STORE YH ON STACK TOP 

00A7 9BD8E1 


193 

FSUB 

ST, ST ( 1 ) 

? CYH-YL) ON STACK TOP 

OOAA 9BDECA 


19? 

FMULP 

ST(2), ST 

? ( ( XV-XL ) / ( XH-XL ) ) ( YH-YL ) ON STACK BOT 

OOAO 9BDEC1 


200 

FADD 


;YV ON STACK TOP (ANSWER) 

OOBO 37DA 


201 

XCHG 

BX, DX 


00B2 EB0B90 


202 

JKP 

SAVEJEILXPTR 




203 




00B5 


204 

OUT-OF-RANGE: 



00B5 9BDDB0 


205 

FST 

ST < 0 ) 

! DISCARD TOP CP STACK 

00B3 8703 


206 

XCHG 

DX, BX 


OOBA 9BD901 


207 

FLD 

DWORD PTR EBX + DI3 

? SAVE FOUNDRY Y VALLE 

OOBO 8703 


203 

XCHG 

DX, BX 




20? 




OOBF 


210 

SAVEJBUFTR: 



OOBF D1EF 


211 

SHR 

DI, 1 

5DI = DI / 4 

0X1 01 EF 


212 

SHR 

DI, 1 


OOC3 83EB03 


213 

SUB 

BX, 8 

;BX = XPTR LOCATION 

00C6 393F 


214 

MOV 

WORD PTR EBX3, DI 

;SAVE NEW XPTR 



215 





12 



3086/37/88/186 MACRO ASSEMBLER RFUN1S 


08/3-0/84 PAGE 5 


LOC OBJ 

LINE 

SOURCE 


00C8 

216 

EPILOGUE: 


00C8 5D 

217 

POP 

BP ; RESTORE REGISTERS 

00C9 C20100 

218 

RET 

1 


21? 




220 

NEARCODE 

ENDP 


221 



— 

222 

CODE ENDS 



223 




224 



225 



226 




227 

END 



ASSEMBLY COMPLETE, NO ERRORS FOUND 
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8036/87/88/186 MACRO ASSEMBLER RFUN1L 


03/30/84 PAGE 1 


SERIES-III 8086/37/33/136 MACRO ASSEMBLER V2.0 ASSEMBLY OF MODULE RFUN1L 
OBJECT MODULE PLACED IN : FI : RFUN1L. OBJ 
ASSEMBLER INVOKED BY: ASM36.36 :F1:RFUN1L.SRC 


LOG OBJ LINE SOURCE 


1 

2 

3 

4 

5 

6 
7 
3 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 




* ROUTINE NAME: RFUN1L * 

* DATE: AUG. 25, 1933 * 

* PROGRAMMER'S NAME: MICHAEL MACK IN * 

* PURPOSE: TO INTEPOLATE A SIMPLE * 

* FUNCTION FOR LARGE-MODEL * 

* PROGRAMS. * 




THIS SUBROUTINE OUTPUTS THE ESTIMATED VALUE OF A FUNCTION YV = F(XV) 
BY INTERPOLATING FROM KNOWN SOLUTION POINTS, USING THE RELATION 
YV = ( ( X V-XL ) / ( XH_XL ) ) { YL-YL ) + YL WHERE 
YV = ANSWER TO BE COMPUTED 
xy = INPUT VALLE 
XH = KNOWN HIGH X POINT 
XL = KNOWN LOW X POINT 
YH = KNOWN HIGH Y POINT 
YL = KNOWN LOW Y POINT 


21 

22 

23 

24 


INPUT REQUIREMENTS: 

1) THE FOLLOWING ADDRESSES WILL BE PASSED TO THIS SUBROUTINE: 
A. ON TOP OF THE 3036 STACK: 



25 

THE ADDRESS 

(OFFSET THEN SEGMENT) OF A DATA AREA WITH 


26 

THE 

FOLLOWING FORMAT 



27 





— 

23 

DSTRUC 

STRUC 




29 





0000 

30 

XPTR 

DW 

7 

JX ARRAY INDEX 

0002 

31 

ZPTR 

DU 

o 

;Y ARRAY INDEX 

0004 

32 

NXPTS 

DW 

7 

5 NO. OF ELEMENTS IN X ARRAY 

0006 

33 

NZPTS 

DW 

o 

5 NO. OF ELEMENTS IN Y ARRAY 


34 i 

; X ARRAY 

DD 

NXPTS DUP<?) 

?X ARRAY 


35 

5 YARRAY 

DD 

NYPT3 DUPC?) 

;Y ARRAY 


36 





— 

37 

DSTRUC 

ENDS 




33 

* 





39 

i B. 2ND FROM THE 

TOP OF THE 3086 STACK: 


40 

41 


THE ADDRESS OF THE MEMORY LOCATION WHICH HOLDS THE X 
INPUT VALUE 


42 

43 

44 

45 

46 

47 
43 

49 

50 


2) IT IS REQUIRED THAT THE 8037 CHIP HAVE 3 EMPTY REGISTERS 
WHEN THIS INTERPOLATION ROUTINE IS CALLED. 

OUTPUT EFFECTS: 

1) REGISTERS DESTROYED: AX, BX, CX, DX. DI, SI 

2) INTERPOLATED RESULT RETURNED ON TOP OF 8037 STACK 

3) XPTR LOCATION UPDATED TO INDEX X VALLE JUST PRECEDING XV. 
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8086/87/38/186 MACRO ASSEMBLER RFUN1L 


08/30/34 PAGE 2 


LOC OBJ 

LINE 

SOURCE 


51 

; EXAMPLE FORTRAN CALL: 


52 

J 


53 

5 SUBROUTINE RTEST 


54 

5 DIMENSION XARRAY<7), YARRAYI7) 


55 

! INTEGER XINBEX. YINDEX 


56 

! COMMON /STUFF/XINDEX, Z INDEX, NXPTS, NZPTS, XARRAYi YARRAY 


57 

! DATA X INDEX . Z INDEX > NXPTSi NZPTS/0. 0. 7 , 7 / 


53 

i DATA XARRAY/O.O. l.Oi 2.0.3. Oi 4. 0i5. 0.6.0/ 


59 

» DATA YARRAY/O.O, 1,0.4.0.9.0.16.0.25.0.36.0/ 


60 

! X=1.5 


61 

> Y=RFUN1L ( X * X INDEX ) 


62 

! PRINT *, Y 


63 

; END 


64 



65 



66 



67 



68 

NAME RFUN1L 


69 

CGROUP GROUP CODE 


70 

DGROUP GROUP DATA 


71 

ASSUME CS: COROUP. DS: [GROUP, ES.’DGROUP 


72 

PUBLIC RFUN1L 


73 



74 

;#M*m**M»M*****««WW*»MM***J««»H»**i#*J**4HHHHHHH«»**M»»# 


75 



76 

; STACK STRUCTURE 


77 


— 

73 

PARAMS STRUC 


79 


0000 

80 

OLD.ES DH ? ; SAVED ES REGISTER 

0002 

31 

OLD.BP DH ? ; SAVED BP REGISTER 

0004 

32 

OLB.DS lW ? ; SAVED DS REGISTER 

0006 

83 

RETURN DD ? ; RETURN ADDRESS (OFFSET ON TOP, THEN SEGMENT) 

OOOA 

84 

MDDR DD ? ; ADDRESS CF DATA STRUCTURE (OFFSET, THEN SEGMENT) 

OOOE 

35 

W CD ? ; ADDRESS OF XV (OFFSET, THEN SEGMENT) 


36 


— 

37 

PARAMS ENDS 


38 



39 



90 



91 


— 

92 

DATA SEGMENT PUBLIC 'DATA' 


93 


0000 ???? 

94 

STHSAVE DH ? 


95 


— 

96 

DATA ENDS 


97 



98 

iJWHHIHmTOHmWHWHHHiJHHWHHIHHIHHHIHHH 


99 

100 


— 

101 

CODE SEGMENT PUBLIC 'CODE' 


102 


0000 

103 

FARCODE PROC FAR 


104 



105 



15 



8036/87/83/ 186 MACRO ASSEMBLER RFUN1L 


08/30/84 PAGE 3 


LQC OBJ 


LINE 

SOURCE 






106 





0000 


107 

RFUN1L: 




0000 IE 


108 

PUSH 

DS 

iSAVE REGISTERS 

0001 55 


109 

PUSH 

BP 



0002 06 


110 

PUSH 

ES 



0003 8BEC 


111 

NOV 

BP, SP 



0005 B3 

R 

112 

MOV 

AX, DATA 



0003 3EC0 


113 

MOV 

ES, AX 



COCA C55E0E 


114 

LDS 

BX, CBP3.XV 

;GET XV ADDRESS AND SEGMENT 

GOOD 9BD907 


115 

FLB 

DWORD PTR EBX3 

if 

5 UT XV ON TOP OF 8087 STACK 

0*010 C55E0A 


116 

LDS 

BX, EBP] . DADDR 

;( 

GET START ADDRESS OF DATA STRUCTURE 

0013 SB3F 


117 

MOV 

DI, EBXLXPTR 

iDI = ELEMENT INDEX 

0015 D1E7 


lie 

SHL 

DI, 1 



0017 B1E7 


ii? 

SHL 

DI, 1 

!DI = BYTE INDEX = ELEMENT INDEX * 4 

001? 8B7704 


120 

MOV 

SI, CBX3.NXPTS 

iSI = NXPTS 

001C B1E6 


121 

SHL 

SI, 1 



001E D1E6 


122 

SHL 

SI, 1 

;SI = NO. OF BYTES=NXPTS * 4 . 

0020 3X303 


123 

ADD 

BX, 3 

? 

BX = START OF X ARRAY 

0023 8BD3 


124 

MOV 

DX, BX 



0025 03D6 


125 

ADD 

DX, SI 

, 

DX = START OF Y ARRAY 

0027 9BB311 


126 +2 

FCGM 

DWORD PTR EBX + DI3 

? 

COMPARE ARRAY ELEMENT WITH XV 

002A 9B26DD3E0000 

R 

127 +1 

FSTSW 

ESsSTWSAVE 



0030 9B 


128 +1 

FWAIT 




0031 26A10000 

R 

129 +1 

MOV 

AX, ESsSTWSAVE 



0035 250041 


130 +1 

AND 

AX, 4100H 

, 

MASK-IN COMPARISON BITS 

0033 3D0001 


131 +2 

CMP 

AX, 0100H 



003B 7434 


132 +2 

JE 

LEFT.Cf .START 


; EXIT IF XV < ARRAY ELEMENT 



133 +2 







134 +1 







135 




IF X > DATAtORIGINAL INDEX) 

003D 


136 

RIGHT JjF.START: 



»** THEN BEGIN 

003D 83C704 


137 

ADD 

DI, 4 


INDEX = INDEX + 4 

0040 3BFE 


133 

RSIs CMP 

DI, SI 


WHILE (INDEX <= NO. OF BYTES) 

0042 7420 


139 

JE 

LI 





140 




Hf [S3 BEGIN «* 

0044 9BB311 


141 +2 

FCGM 

DWORD PTR [BX + DI3 


COMPARE ARRAY ELEMENT WITH XV 

0047 9B26DD3EOGOO 

R 

142 +1 

FSTSW 

ESsSTWSAVE 



004D 9B 


143 +1 

FWAIT 




004 E 26A10000 

R 

144 +1 

MOV 

AX, ESsSTWSAVE 



0052 250041 


145 +1 

AND 

AX, 4100H 

' 

MASK-IN COMPARISON BITS 

0055 300001 


146 +2 

CMP 

AX, 0100H 



0053 740A 


147 +2 

JE 

LI 

1 

EXIT IF XV < ARRAY ELEMENT 

005A 3B0040 


148 +2 

CMP 

AX, 4000H 



005B 7469 


14? +2 

JE 

BREAKPOINT 


! EXIT IF XV = ARRAY ELEMENT 



150 +1 





005F 33C704 


151 

ADD 

DI, 4 


INDEX = INDEX + 4 

0062 EBDC 


152 

JMP 

RSI 


*** END (WHILE) *** 

0064 83EE04 


153 

Lis SUB 

SI, 4 



0067 S3EF04 


154 

SUB 

DI. 4 



006A 3BFE 


155 

CMP 

DI, SI 


IF INDEX > NO. Cf BYTES 

006 C 745A 


156 

JE 

OUT-OF-RANGE 


THEN 0LOF.RANOE 

006E EB3490 


157 

JMP 

FGUND_LOW_ INDEX 


ELSE FOUNB-LOILIHDEX 



153 




*** END (IF) *** 



15? 




IF X < DATAtORIGINAL INDEX) 

0071 


160 

LEFT.OF.START: 



*»* THEN BEGIN *»* 


16 



8086/87/83/186 MACRO ASSEMBLER RFUN1L 


03/30/84 FACE 4 


LOC OBJ 


LINE 

SOURCE 



0071 83EF04 


161 

SUB 

Bit 4 

! INDEX = INDEX - 4 

0074 83FF00 


162 

LSI; CMP 

DI, 0 

! WHILE (INDEX > 0) 

0077 7C20 


163 

JL 

Ml 

* 



164 



i »** DO BEGIN *** 

007? 9BD811 


165 +2 

fcom 

DWORD PTR EBX + DII 

i COMPARE ARRAY ELEMENT WITH XV 

007C 9B260D3E0000 

R 

166 +1 

FSTSW 

ES.'STHSAVE 


0082 9B 


167 +1 

FIJAIT 



0033 26A10000 

R 

163 +1 

MOV 

AX, ESiSTHSAVE 


0087 250041 


16? +1 

ANB 

AX. 4100H 

' MASK-IN COMPARISON BITS 

003A 3DOOOO 


170 +2 

CMP 

AX, OOOOH 

! EXIT IF XV > ARRAY ELEMENT 

OOSD 740A 


171 +2 

JE 

Ml 


003F 3D0040 


172 +2 

CMP 

AX, 4000H 


0092 7434 


173 +2 

174 +1 

JE 

BREAKPOINT 

! EXIT IF XV = ARRAY ELEMENT 

0094 83EF04 


175 

SUB 

DI, 4 

! INDEX = INDEX - 4 

0097 EBDB 


176 

JHP 

LSI 

*** END (WHILE) *♦* 

0099 S3FF00 


177 

Ml: CMP 

DI, 0 

? IF INDEX O 0 

0090 7D06 


173 

JOE 

FOUND J.QH-INDEX 

! THEN FOUND_LQW_ INDEX 

009E 83C704 


17? 

ABB 

DI, 4 

! 

00 A 1 EB2590 


130 

JMP 

OUT-OF-RANGE 

! ELSE OUT-OF_RANGE; 



181 



;#** END (ID) *4* 



132 




00A4 


183 

FOUND_L0W-INDEX: 


00A4 


184 

INTERPOLATE: 



00A4 9BB901 


185 

FLB 

DWORD PTR [BX+DI3 

; STORE XL ON STACK TOP 

00A7 9BBCE9 


186 

FSUB 

ST ( 1 > t ST 

» ( XV— XL ) ON STACK BOTTOM 

OOAA 9BB86904 


187 

FRIER 

DWORD PTR (BX+DI+4) 

? (XH-XL) ON STACK TOP 

OOAE 9BDEF9 


183 

FBIV 


J (XV— XL)/(XH— XL) LEFT ON STACK 

GOBI 87BA 


18? 

XCHG 

BX, DX 


00B3 9BD901 


190 

FLB 

DWORD PTR CBX+DI3 

; STORE YL ON STACK TOP 

00B6 9BB94104 


191 

FLB 

DWORD PTR CBX+DI+43 

? STORE YH ON STACK TOP 

OOBA 9BB3E1 


192 

FSUB 

ST, ST'l) 

; (YH-YL) ON STACK TOP 

OOBB 9BBECA 


193 

FMULP 

ST(2), ST 

j ( (XV— XL)/ (XH-XL) ) (YH-YL)QN STACK EOT 

COCO 9BDEC1 


194 

FABB 


,-YV ON STACK TOPf ANSWER) 

00C3 87BA 


195 

XCHG 

BX, BX 


0005 EB0B90 


196 

JMP 

SAUEJBLXPTR 




197 




0008 


193 

BREAKPOINT: 



OOC8 


199 

OUT.OF-RANGE: 



0003 9BDBD8 


200 

FSTP 

ST(O) 

5 DISCARD TOP OF STACK 

OOCB S7B3 


201 

XCHG 

BX, BX 


OOCD 9BD901 


202 

FLB 

DWORD PTR EBX + DI3 

5SAVE BOUNDRY Y VALLE 

OOBO 87B3 


203 

XCHG 

DX, BX 




204 




00B2 


205 

SAVE_NEH_XPTR: 



OOB2 B1EF 


206 

SHR 

DI, 1 

!DI = DI / 4 

00B4 B1EF 


207 

SHR 

DI, 1 


00B6 83EB08 


208 

sub 

BX, 8 

,'BX = XPTR LOCATION 

COD? 893F 


20? 

MOV 

WORD PTR [BX], DI 

.'SAVE NEW XPTR 



210 




OOBB 


211 

EPILOGUE: 



OOBB 07 


212 

POP 

ES 


OOBC 5B 


213 

POP 

BP 

! RESTORE REGISTERS 

OOBB IF 


214 

POP 

BS 


OOBE CA0300 


215 

RET 

8 



17 



8036/87/88/186 MACRO ASSEMBLER RFUN1L 


08/30/84 PAGE 5 


LOC OBJ LINE SOURCE 

216 

217 FARCODE ENBP 

218 

219 CODE ENDS 

220 

223 

224 END 
ASSEMBLY COMPLETE, NO ERRORS FOUND 
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8086/87/83/186 MACRO ASSEMBLER RFUN2S 


03/30/84 PAGE 1 


SERIES-III 3036/87/83/186 MACRO ASSEMBLER V2.0 ASSEMBLY OF MODULE RFUN2S 
OBJECT MODULE PLACED IN :Fl:RFUN2S.0BJ 
ASSEMBLER IN'/OKED BY: ASM36.86 :Fl:RFUN2S.SRC 


LOC OBJ LINE SOURCE 


mMm«m#*w****»#**###*#»*m****** 


* ROUTINE NAME: RFUN2S * 

* DATE: AUG. 27, 1983 * 

* PROGRAMMER'S NAME: MICHAEL MACKIN * 

* PURPOSE: TO INTERPOLATE A SIMPLE * 

* FUNCTION FOR SMALL-MODEL » 

* PROGRAMS * 


THIS SUBROUTINE OUTPUTS THE ESTIMATED VALUE OF A FUNCTION YV = F(XV.ZV) 
BY INTERPOLATING FROM KNOWN SOLUTION POINTS, USING THE RELATIONS 
YL = ( ( XV-XL ) / ( XH— XL ) ) ( YB-YA) + YA 
YH = ( ( XV-XL ) / ( XH-XL) ) ( YD-YC) + YC 
YV = ( ( ZV-ZL ) / ( ZH.ZL ) ) ( YH-YL ) + YL 

INFUT REQUIREMENTS: 

1) THE FOLLOWING ADDRESSES WILL BE PASSED TO THIS SUBROUTINE: 

A. ON TOP OF THE 8036 STACK: 

THE ADDRESS (OFFSET) OF A DATA AREA WITH THE FOLLOWING 
FORMAT 

DSTRUC STRUG 
XPTR DvJ ? 


1 

2 

3 

4 

5 

6 

7 

8 
Q 

10 

11 

12 

13 

14 

15 

16 
17 
13 
1? 
20 



21 

— 

22 


23 

0000 

24 

0002 

25 

0004 

26 

0006 

27 


ZO 

2? 

30 


31 

oo 


•Ji. 

33 


34 


35 


36 


37 


33 


3? 


ZPTR 

DM 

9 

NXPTS 

DU 

9 

NZPTS 

DM 

? 

ZARRAY 

DD 

NZPTS DU?(?) 

XARRAY 

DD 

NXPTS DUP(?) 

YARRAY 

DD 

NYPTS DUP(?) 

DSTRUC 

ENDS 



B. ON THE 8037 STACK: 
ON TOP ZV, THEN XV 


JX ARRAY INDEX 

;z ARRAY INDEX 

SNQ. OF ELEMENTS IN X ARRAY 

5 NO. OF ELEMENTS IN Z ARRAY 

5Z ARRAY (255 ELEMENTS MAXIMUM) 

?X ARRAY 

5Y ARRAY 


2) IT IS REQUIRED THAT THE 8087 CHIP HAVE 5 EMPTY REGISTERS 
UHEN THIS INTERPOLATION ROUTINE IS CALLED. 


40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


3) NOTE: NO THO ADJACENT ELEMENTS OF XARRAY GR YARRAY MAY HAVE 
IDENTICAL ELEMENTS. THIS HILL RESULT IN DIVISION BY ZERO. 

OUTPUT EFFECTS: 

1) REGISTERS DESTROYED: AX, BX, CX, DX, DI, SI 

2) INTERPOLATED RESULT RETURNED ON TOP OF 3087 STACK 

3 ) XPTR LOCATION UPDATED TO INDEX X VALUE JUST PRECEDING XV. 
ZPTR LOCATION UPDATED TO INDEX Z VALUE JUST PRECEDING ZV. 
(TO SPEED FUTURE SEARCHES FOR X AND Z) 

EXAMPLE PL/M CALL: 


19 



3036/37/83/186 
LOG OBJ 


0000 


MACRO ASSEMBLER RFUN2S 


03/30/34 PAGE 2 


LINE SOURCE 
51 


52 5 


‘SMALL 


53 ; 




54 ; 


TSTFGM: 

DO; 

55 




56 ; 


/* THIS IS A PROGRAM TO TEST ROUTINE RFUN2S */ 

57 ; 




58 



DECLARE DSTRUC STRUCTURE ( 

5? 



XINDEX INTEGER* 2INDEX INTEGER* NXPTS INTEGER* NZPTS INTEGER 

60 



ZARRAYC3) REAL* X ARRAY! 5) REAL* YARRAY(15) REAL) PUBLIC 

61 



INITIAL ( 

62 



/* XINDEX */ 0* 

63 



/* Z INDEX */ 0* 

64 



/* NXPTS */ 5* 

65 



/* NZPTS */ 3* 

66 



/* ZARRAY */ 0.0* 10.0, 20.0, 

67 



/* XARRAY */ 0.0* 2.0* 4.0* 6.0* 3.0* 

68 



/* YARRAY */ 0.0* 2.0* 4.0, 6.0* 3.0, 

6? 



10.0* 12.0, 14.0* 16.0* 13.0* 

70 



20.0, 22.0* 24.0* 26.0* 23.0 )J 

71 




72 



RFUN2S: PROCEDURE ( X , Z , STRUC. ADDR ) REAL EXTERNAL; 

73 



DECLARE X REAL; 

74 



DECLARE Z REAL; 

75 



DECLARE STRUCLADDR POINTER; 

76 



END RFUN2S; 

77 




78 



FUNPGM: PROCEDURE PUBLIC? 

7? 



DECLARE I INTEGER; 

80 I 



DECLARE (X* Z,Y) REAL? 

31 



X = 1.5? 

82 



o 

CN 

II 

r-M 

33 



Y = RFUN2S ( X , Z * 6DSTRUC > * 

84 



CALL PRINT-ANSWER; 

35 



END; 

86 



END FUNPGM; 

87 




83 


END TSTPGM? 

3? 




80 


?1 




92 




83 


NAME 

RFUN2S 

84 

CGROUP 

GROUP 

CODE 

85 

DGROUP 

GROUP 

DATA 

86 


ASSUME 

CS: CGROUP* DS: DGROUP 

87 


PUBLIC 

RFUN2S 

83 




88 


100 




101 

» 

STACK STRUCTURE 

102 




103 

PARAKS 

STRUC 


104 




105 

OLD-BP 

DW 

? i SAVED BP REGISTER 


20 



8086/37/83/186 MACRO ASSEMBLER RFUN2S 


03/30/84 PAGE 3 


LOG OBJ 


LINE 

SOURCE 


0002 


106 

RETURN BW 

? ; RETURN ADDRESS 

0004 


107 

DADDR DM 

? ! ADDRESS OF DATA STRUCTURE 



103 



— 


10? 

PARAMS ENDS 




110 





111 




112 




113 



— 


114 

DATA SEGMENT PUBLIC / [iATA / 



115 



0000 ???? 


116 

STHSAVE DM 

9 



117 



— 


113 

DATA ENDS 




11? 





120 




121 




122 



— 


123 

CODE SEGMENT PUBLIC "CODE" 



124 



0000 


125 

NEARCOBE 

PROC NEAR 



126 





127 





123 



0000 


12? 

RFUN2S* 




130 


; SAVE REGISTERS 

0000 55 


131 

PUSH 

BP 

0001 SBEC 


132 

MOV 

BP, SP 

0003 3B5E04 


133 

MOV 

BX, [DPI. DADDR 5 GET START ADDRESS Cf DATA STRUCTURE 

0006 SBB3 


134 

MOV 

DX, BX 5 SAVE BX 

0003 3B0F 


135 

MOV 

CX, IBXLXPTR SCX = XPTR 

OOOA SB6F04 


136 

MOV 

BP, CBXJ.NXPTS 5 BP = NXPTS 

GOOD 8B7F02 


137 

MOV 

DIi IBXLZPTR ?DI = ZPTR 

0010 D1E7 


133 

SHL 

DI, 1 

0012 B1E7 


13? 

SHL 

DI. 1 5DI = BYTE INDEX = ELEMENT INDEX * 4 

0014 SB7706 


140 

MOV 

SI, [BX]. NZPTS ;SI = NZPTS 

0017 D1E6 


141 

SHL 

SI, 1 

001? D1E6 


142 

SHL 

SI, 1 ;SI = NO. OF BYTES = NZPTS * 4 



143 



001 B 


144 

GET.Z.INDEX: 


00 IB 83C308 


145 

ADD 

BX, 8 ;BX = START OF Z ARRAY 

001 E 9BB311 


146 +2 

FCOM 

DWORD PTR [BX + DI] 5 COMPARE ARRAY ELEMENT WITH ZV 

0021 9BDB3E0000 

R 

147 +1 

FSTSW 

STHSAVE 

0026 ?B 


143 +1 

FNAIT 


0027 ftlOOOO 

R 

14? +1 

MOV 

AX, STHSAVE 

002A 250041 


150 +1 

AND 

AX> 4100H ; MASK-IN COMPARISON BITS 

«X?2B 3B0001 


151 +2 

CMP 

AX, 0100H 

0030 7432 


152 +2 

JE 

LEFT.OF.Z.START ; EXIT IF ZV < A 




MENT 




153 +2 





154 


! IF Z > DATAIORIGINAL INDEX) 

0032 


155 

RIGHT.0F_2_START: ;##» THEN BEGIN *** 

0032 83C704 


156 

ADD 

DI, 4 i INDEX = INDEX + 4 

0035 3BFE 


157 

RZSl: CMP 

DI, SI ; WHILE (INDEX <C= NO. OF BYTES) 

0037 741? 


158 

JE 

RZS2 ; 



15? 


! ♦** DO BEGIN *** 


21 



S0S6/S7/SS/ 1 86 MACRO ASSEMBLER RFUN2S 


08/30/84 PAGE 4 


LOC OBJ 


LINE 

SOURCE 



003? 9BD311 


160 +2 


FCON 

DWORD PTR EBX + Dll 

003C 9PDD3EOOOO 

R 

161 +1 


FST3W 

STWSAVE 

0041 ?B 


162 +1 


FIJAIT 


0042 A 10000 

R 

163 +1 


MOV 

AX 7 STWSAVE 

0045 250041 


164 +1 


AND 

AX > 4100H 

0048 3D0001 


165 +2 


CMP 

AX, 010CH 

004B 7405 


166 +2 


JE 

RZS2 



167 +2 




004D 83C704 


168 


ADD 

DI. 4 

0050 EBE3 


16? 


JMP 

RZS1 

0052 8BC6 


170 

RZS2: 

MOV 

AX 7 SI 

0054 2DM00 


171 


SUB 

AX, 4 

0057 83EF04 


172 


SUB 

DI, 4 

005A 3EF8 


173 


CMP 

DI, AX 

005C 7540 


174 


JNE 

FOUNElLOILZ. INDEX 

005E 9BDDD3 


175 


FSTP 

ST ( 0 ) 

0061 EB4B90 


176 


JMP 

GETJLIKDEX 



177 






178 




0064 


17? 

LEFOF, 

-Z-START 

: 

0C>64 83EF04 


180 


SUB 

DI, 4 

0067 83FFOO 


131 

LZSi: 

CMP 

DI, 0 

006A 7C1E 


182 


JL 

LZS2 



133 




006C 9BD311 


134 +2 


FCOM 

DWORD PTR CBX + DI3 

006F 9BDD3E0000 

R 

135 +1 


FSTSW 

STWSAVE 

0074 9B 


186 +1 


FWAIT 


0075 A 10000 

R 

187 +1 


MOV 

AX, STWSAVE 

0078 250041 


183 +1 


AND 

AX, 4100H 

007B 3DOOOO 


IS? +2 


CMP 

AX, OOOOH 

007E 740A 


190 +2 


JE 

LZS2 

0080 3D0040 


191 +2 


CMP 

AX, 4000H 

0083 7405 


192 +2 


JE 

LZS2 

0035 33EF04 


193 


SUB 

DI, 4 

0083 EBDD 


194 


JMP 

LZSI 

008A 83FFOO 


195 

LZS'2: 

CMP 

DI, 0 

0080 7D0F 


196 


JGE 

FQUNELLOW-ZJNDEX 

008F 33C704 


197 


ADD 

DI, 4 

0092 9BDDD3 


198 


FSTP 

ST(O) 

0095 9BD9EE 


19? 


FLDZ 


0098 9BD9C? 


200 


FXCH 

ST ( 1 ) 

009B EB1190 


201 


JMP 

GET.X.INDEX 



202 






203 




009E 


204 

FOUNB-LOW-ZJNDEX: 

009E 9BD901 


205 


FLD 

DWORD PTR CBX+DI) 

00 A 1 9BDCE? 


206 


FSUB 

ST ( 1 ) , ST 

00A4 9BD36904 


207 


FSUBR 

DWORD PTR (BX+DI+4] 

WAS 9BDEF9 


208 


FDIV 


OOAB 9BD9C9 


20? 


FXCH 

ST(1> 



210 




OOAE 


211 

GETJLINDEX: 


OOAE 03DE 


212 


ADD 

BX, SI 

00 BO 87F? 


213 


XCHG 

DI, CX 

00B2 D1E7 


214 


SHL 

DI, 1 


COMPARE ARRAY ELEMENT WITH ZV 


MASK-IN COMPARISON BITS 

EXIT IF ZV < ARRAY ELEMENT 

INDEX = INDEX + 4 
*** END (WHILE) *** 


IF INDEX O NO. OF BYTES 
THEN F0UND.L01L INDEX 
ELSE 0UT_0F-RANGE 

*#* END (IF) 

IF Z < DATA (ORIGINAL INDEX) 

*** THEN BEGIN *** 

INDEX = INDEX - 4 
WHILE (INDEX > 0) 

##* DO BEGIN *** 

COMPARE ARRAY ELEMENT WITH ZV 


MASK-IN COMPARISON BITS 
EXIT IF ZV > ARRAY ELEMENT 


EXIT IF ZV = ARRAY ELEMENT 
INDEX = INDEX - 4 
*** END (WHILE) *** 

IF INDEX O 0 

THEN FOUNELLOW.INDEX 
ELSE OUT-GF.RANGE 

REPLACE TOP WITH 0 
PUT XV ON TOP OF 8087 ST. 

*** END{IF) **# 


; STORE ZL ON STACK TOP 
5 (ZV-ZL) 2ND FROM TOP 
? ( ZH— ZL) ON STACK TCP 
; (ZV-ZL) /(ZH-ZL) LEFT ON STACK 
;PUT XV ON TOP OF 8087 STACK 


?BX = START OF X ARRAY 

;DI = XPTR : CX = LOHLZJNDEX 
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5 


LOC OBJ 


LINE 

SOURCE 



00E4 B1E7 


215 


SHL 

DI, 1 

00B6 37F5 


216 


XCh‘G 

SI, BP 

00B8 D1E6 


217 


SHL 

SI, 1 

OOBA D1E6 


213 


SHL 

SI, 1 

ODBC 9BD31 1 


219 +2 


FCOM 

CHORD PTR CBX + DI] 

OOBF 9BDD3EQOOO 

R 

220 +1 


FSTSW 

STHSAVE 

00C4 9B 


221 +1 


FHAIT 


ooc5 mm 

R 

222 +1 


MOV 

AX, STUSAVE 

OOCS 250041 


223 +1 


AND 

AX, 4100H 

OOCB 3B0001 


224 +2 


CMP 

AX, 0100H 

OOCE 7432 


225 +2 


JE 

LEFT.OF-X.START 




WENT 





226 +2 






227 




OODO 


223 

RIGHT.OF.X.START: 

OODO 83C704 


229 


ADD 

DI, 4 

00D3 3BFE 


230 

RXSl: 

CMP 

DI, SI 

00D5 741? 


231 


JE 

RXS2 



232 




OOD7 9BD811 


233 +2 


FCOM 

DWORD PTR [BX + DI3 

OODA 9BDB3E0000 

R 

234 +1 


FSTSW 

STUSAVE 

CODE 9B 


235 +1 


FHAIT 


OOEO A 10000 

R 

236 +1 


MOV 

AX, STUSAVE 

COE 3 250041 


237 +1 


AND 

AX, 4100H 

00E6 3B0001 


238 +2 


CMP 

AX, 0100H 

00E9 7405 


23? +2 


JE 

RXS2 



240 +2 




OOEB 83C704 


241 


ADD 

DI, 4 

00 EE EPE3 


242 


JMP 

RXSl 

OOFO SBC6 


243 

RXS2: 

MOV 

AX, SI 

00F2 2B0400 


244 


SUB 

AX, 4 

OOFS 83EF04 


245 


SUB 

DI, 4 

OOFS 3BF8 


246 


CMP 

DI, AX 

OOFA 753D 


247 


JNE 

FGUND_LQO_ INDEX 

OOFC 9BD3D3 


243 


FS'TP 

ST(O) 

OOFF EB4590 


24? 


JMP 

GET_YL 



250 






251 




0102 


252 

LEFT.GF.X.START: 

0102 83EF04 


253 


SUB 

DI, 4 

0105 83FF00 


254 

LXS'l : 

CMP 

DI, 0 

0108 7C1E 


255 


JL 

LXS2 



256 




OlOA 9BD811 


257 +2 


FCOM 

DWORD PTR CBX + DIJ 

01 OD 9BDB3E0000 

R 

253 +1 


FSTSW 

STHSAVE 

0112 9B 


25? +1 


FWAIT 


0113 A 10000 

R 

260 +1 


MOV 

AX, STHSAVE 

0116 250041 


261 +1 


AND 

AX, 4100H 

Oil? 3B0000 


262 +2 


CMP 

AX, OOOOH 

011C 740A 


263 +2 


JE 

LXS2 

01 IE 3D0040 


264 +2 


CMP 

AX, 4000H 

0121 7405 


265 +2 


JE 

LXS2 

0123 83EF04 


266 


SUB 

DI, 4 

0126 EBDD 


267 


JMP 

LXS1 

0128 83FF00 


268 

LXS2: 

CMP 

DI, 0 


DI = BYTE INDEX = XPTR * 4 
SI = NXPTS 

SI = NO. OF X BYTES = NXPTS * 4 
COMPARE ARRAY ELEMENT WITH XV 


MASK-IN COMPARISON BITS 

? EXIT IF XV < ARRAY ELE 


IF X > DATA (ORIGINAL INDEX) 

»** THEN BEGIN w 
INDEX = INDEX + 4 
WHILE (INDEX <= NO. OF BYTES) 

«* DO BEGIN *#* 

COMPARE ARRAY ELEMENT WITH XV 


MASK- IN COMPARISON BITS 

EXIT IF XV < ARRAY ELEMENT 

INDEX = INDEX + 4 
*** END (WHILE) *** 


IF INDEX O NO. OF BYTES 
THEN FOUND-LOW. INDEX 
ELSE OUT.C'FJtANGE 

*** END (IF) *** 

IF X < DATA (ORIGINAL INDEX) 

*** THEN BEGIN *** 

INDEX = INDEX - 4 
WHILE (INDEX > 0) 

*** DO BEGIN *** 

COMPARE ARRAY ELEMENT WITH XV 


MASK-IN COMPARISON BITS 
EXIT IF XV > ARRAY ELEMENT 


EXIT IF XV = ARRAY ELEMENT 
INDEX = INDEX - 4 
**» END (WHILE) *** 

IF INDEX O 0 
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LOG OBJ 

LINE 

SOURCE 



012B 7D0C 

269 

JGE 

FQUND-LOHJLINBEX 

THEN FOUND-LOW- INDEX 

01 2D 83C704 

270 

ADD 

DI, 4 

ELSE OUT.OF.RANGE 

0130 9BDDD8 

271 

FSTP 

ST(O) 


0133 9BD9EE 

272 

FLDZ 


REPLACE TOP WITH 0 

0136 EB0E90 

273 

JHP 

GET.YL 



274 



*** END (IF) *** 


275 




0139 

276 

FOUNOJ.GOJNDEX: 


0139 9BD901 

277 

FLD 

CHORD PTR [BX+DI3 5 STORE XL ON STACK TOP 

013C 9BDCE9 

273 

FSUB 

ST { 1 ) » ST 

(XV-XL) ON STACK BOTTOM 

013F 9BD36904 

279 

FSUBR 

DWORD PTR [BX+DI+4] 

(XH-XL) ON STACK TOP 

0143 9BDEF9 

230 

FDIV 


{ XV-XL ) / ( Xf-f-XL ) LEFT ON STACK 


281 




0146 

232 

GET- YL ? 



0146 03DE 

283 

ADD 

BX, SI ;BX = START OF YARRAY 

0148 8EC6 

234 

MOV 

AX, SI ;AX = NO. OF X BYTES 

014A D1E3 

285 

3HR 

AX, 1 


014C D1E3 

236 

SHR 

AX, 1 ?AX = NXPTS = NO. OF X BYTES / 4 

014E F6E1 

287 

MUL 

CL ;AX = YA ROW OFFSET 


283 



THE AVOVE IHSTR. LIMITS LENGTH OF ZARRAY 


28? 



TO 255 ELEMENTS 

0150 03D3 

290 

ADD 

BX, AX 

BX = YA ROW PTR = NXPTS * ZPTR 

0152 03DF 

291 

ADD 

BX, DI 

ex = YA POINTER 

0154 9BD907 

292 

FLD 

CHORD PTR [BX3 

STORE YA ON ST TOP 

0157 SBC? 

293 

MOV 

AX, DI 


0159 050400 

294 

ADD 

AX, 4 


015C 3BC6 

295 

CMP 

AX, SI 


015E 7410 

296 

JE 

GET.YH 5EXIT IF XV PAST FOUNDRY 

0160 9ED9C9 

297 

FXCH 

ST ( 1 ) 

(XV-XL)/! XH-XL) ON TOP, THEN YA 

0163 9BD94704 

OQO 

FLD 

DWORD PTR CBX + 4] 

STORE YB ON ST. TOP 

0167 9BDSE2 

29? 

FSUB 

ST, ST (2) 

(YB-YA) ON ST. TOP 

01 6 A 9BDSC9 

300 

FMUL 

ST, ST < 1 ) 

( ( XV-XL) / ( XH-XL) ) (YB-YA) ON ST TOP 

01 6D 9EDEC2 

301 

FADDP 

ST (2) , ST 

YL 2ND FROM TOP OF STACK 


302 




0170 

303 

GET-YH: 



0170 8BC1 

304 

MOV 

AX, CX 5 IF (Z.POINTER TOO BIG) 

0172 050400 

305 

ADD 

AX, 4 ;*** DO BEGIN *** 

0175 3BC5 

306 

CMP 

AX, BP 


0177 750F 

307 

JNE 

GYHl 


0179 SBC7 

308 

MOV 

AX, DI 

IF (X-F’O INTER TOO BIG) 

017B 050400 

30? 

ADD 

AX, 4 


01 7E 3BC6 

310 

CMP 

AX, SI 


0180 742A 

311 

JE 

FOUND.YV 

THEN EXIT 

0182 9BDDD3 

312 

FSTP 

ST 

ELSE DISCARD JUNK 

0185 EB2590 

313 

JM P 

FOUND.YV 

AND THEN EXIT; 


314 



*** END {IF} *♦* 

0183 03DE 

315 

GYHl: ADD 

BX, SI 

BX = YC POINTER 

018A 9BD907 

316 

FLD 

DWORD PTR [BX] 

STORE YC ON STACK TOP 

018D 8BC7 

317 

MOV 

AX, DI 


01 8F 050400 

318 

ADD 

AX, 4 


0192 3BC6 

31? 

CMP 

AX, SI 


0194 740D 

320 

JE 

GET.YV 

EXIT IF XV PAST BOUNDRY 

0196 9BD94704 

321 

FLD 

DWORD PTR [BX + 43 

STORE YD ON STACK TOP 

019A 9BDSE1 

322 

FSUB 

ST, ST < 1 ) 

;YD-YC ON ST TOP 

019D 9BDECA 

323 

FKULP 

ST ( 2 ) * ST 

( ( XV-XL ) / ( XH-XL ) ) ( YD-YC ) 2ND FROM TOP 
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LOC OBJ 

LINE 

SOURCE 



01 AO 9BBEC1 

324 

FADD 


;YH ON TCP Cf 8037 


325 




01A3 

326 

GET.YV: 


lYH TOP, THEN YL, THEN SLOPE 

01 A3 9BDSE1 

327 

FSUB 

ST, ST( i ) 

?YH-YL ON ST. TOP 

01A6 9BBECA 

323 

FMULP 

ST(2), ST 

5 ( ( Z V-ZL ) / ( ZH-ZL ) ) ( YH-YL ) 2ND FR( 

01A? 9BDEC1 

32? 

FADD 


;vv m st. tcp 


330 




01 AC 

331 

FOUND.W: 



01AC 

332 

SAVE.NEW_POINTERS: 


01 AC SBDA 

333 

MOV 

ex, dx 


01AE B1EF 

334 

SHR 

DI, 1 


01 BO DIEF 

335 

SHR 

DI, 1 


01 B2 893F 

336 

MOV 

EBX3.XPTR, DI 

' SAVE XPTR 

01B4 DIE? 

337 

SHR 

CX, 1 


01B6 DIE? 

333 

SHR 

CX, 1 


01B8 894F02 

33? 

MOV 

EBXI.ZPTR, CX 

I SAVE ZPTR 


340 




01 BB 

341 

EPILOGUE: 



01 BB 5D 

342 

POP 

BP 

5 RESTORE REGISTERS 

01 PC C20100 

343 

RET 

1 



344 





345 

NEARCODE 

ENDP 



346 




— 

347 

CODE ENDS 




34S 





34? 

)******MMM********MM**M*MM*M**M**MM*MmmM************ 


350 



351 





352 

END 




ASSEMBLY COMPLETE, NO ERRORS FOUND 
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SERIES— III 8036/37/83/186 MACRO ASSEMBLER V2.0 ASSEMBLY OF MODULE RFUN2L 
OB.JECT MODULE PLACED IN :F1>RFUN2L.0BJ 
ASSEMBLER INVOKED BY: ASM36.36 SFURFUM2L.SRC 


LOC OBJ LINE SOURCE 


**M**mMMM*#**#»**»**»**##**»******* 


* ROUTINE NAME: RFUN2L * 

* DATE: SEP. 1, 1983 * 

* PROGRAMMER'S NAME: MICHAEL MACKIN * 

* PURPOSE: TO INTERPOLATE A SIMPLE * 

* FUNCTION FOR LARGE-MODEL * 

* PROGRAMS * 


THIS SUBROUTINE OUTPUTS THE ESTIMATED VALUE OF A FUNCTION YV = F(XV.ZV) 
BY INTERPOLATING FROM KNOWN SOLUTION POINTS. USING THE RELATIONS 
YL = ( ( X V-XL ) / ( XH-XL ) ) ( YB-Y A ) + YA 
YH = ( ( X V-XL )/( XH-XL )) ( YD-YC ) + YC 
YV = ( ( ZV-ZL ) / ( ZH.ZL ) ) ( YH-YL ) + YL 

INPUT REQUIREMENTS: 

1) THE FOLLOWING ADDRESSES HILL BE PASSED TO THIS SUBROUTINE: 

A. ON TOP OF THE 8086 STACK: 

THE ADDRESS (OFFSET) OF A DATA AREA WITH TIE FOLLOWING 
FORMAT 

DSTRUC STRUC 


1 

2 

3 

4 

5 

6 
7 
3 

9 

10 

11 

12 



13 


14 


15 


16 


17 


13 


19 


20 


21 

CO 


Li 

23 

0000 

24 

0002 

25 

0004 

26 

0006 

27 


23 


29 


0 J 

31 

— 

32 


33 


34 


35 


36 


37 


33 


39 


40 


41 


42 


43 


44 


45 


46 


47 


48 


49 


50 


XPTR 

Dl-J 

7 

ZPTR 

m 

7 

NXPTS 

m 

0 

NZPTS 

m 

? 

ZARRAY 

DD 

NZPTS BUP(?) 

XARRAY 

DD 

NXPTS DUP ( ? > 

YARRAY 

DD 

NYPTS BUP(?) 

DSTRUC 

ENDS 


B. NEXT 

ON THE 8086 STACK: 


;X ARRAY INDEX 
M ARRAY INDEX 
5 NO. OF ELEMENTS IN X ARRAY 
5W0. OF ELEMENTS IN 2 ARRAY 
;Z ARRAY (255 ELEMENTS MAXIMUM) 
JX ARRAY 


JY ARRAY 


THE ADDRESS OF THE MEMORY LOCATION WHICH HOLDS THE X 
INPUT VALUE. 

2) IT IS REQUIRED THAT THE 8037 CHIP HAVE 5 EMPTY REGISTERS 
WHEN THIS INTERPOLATION ROUTINE IS CALLED. 

3) NOTES NO TWO ADJACENT ELEMENTS OF XARRAY OR YARRAY MAY HAVE 
IDENTICAL ELEMENTS. THIS WILL RESULT IN DIVISION BY ZERO. 

OUTPUT EFFECTS: 

1) REGISTERS DESTROYED: AX, BX, CX, DX, DI, SI 

2) INTERPOLATED RESULT RETURNED ON TOP OF S037 STACK 

3) XPTR LOCATION UPDATED TO INDEX X VALUE JUST PRECEDING XV. 
ZPTR LOCATION UPDATED TO INDEX Z VALUE JUST PRECEDING ZV. 
(TO SPEED FUTURE SEARCHES FOR X AND Z) 
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8086/87/83/186 MACRO ASSEMBLER 


LGC OBJ 

LINE 

SOURCE 


51 


EXAMPLE FORTRAN CALL 


52 




53 

; 

SUBROUTINE TEST 


54 


DIMENSION XARRAY(5), ZARRAY<3), YARRAYU5) 


55 

! 

INTEGER XINDEX, ZINDEX, NXPTS, NZPTS 


56 


COMMON /STUFF/XINDEX, ZINDEX, NXPTS, NZPTS, ZARRAY, X ARRAY, YARRAY 


57 

: 

DATA XINDEX, ZINDEX, NXPTS, NZPTS/0, 0,5, 3/ 


53 


DATA ZARRAY/O. 0,10.0,20.0/ 


59 

5 

DATA XARRAY/0. 0, 2. 0, 4. 0, 6. 0, 8. 0/ 


60 

i 

DATA YARRAY/O. 0, 2. , 4. , 6. , 8. , 10. , 12. , 14. , 16. , 18. ,20. ,22. , 24. , 26, 


61 

i 

X=2.5 


62 

' 

Z=12.0 


63 


Y=RFUN2L(X,Z, XINDEX) 


64 

! 

PRINT *, Y 


65 


END 


66 




67 

i 



63 




69 




70 


NAME RFUN2L 


71 

CGROUP GROUP CODE 


72 

DGROUP GROUP DATA 


73 


ASSUME CSiCGROUPi DS:BGROLP» ES.'DGROUP 


74 


PUBLIC RFUN2L 


75 




76 

1 



77 




78 

J 

STACK STRUCTURE 


7? 



— 

30 

FARMS STRUC 


81 



0000 

82 

OLD.ES DM ? ; SAVED ES REGISTER 

0002 

S3 

OLD.BP DH ? ; SAVED BP REGISTER 

0004 

84 

OLD.DS DM ? * SA'.'ED DS REGISTER 

0006 

85 

RETURN DD ? i RETURN ADDRESS (OFFSET ON TOP, THEN SEGTENT) 

OOOA 

86 

DADDR DD ? ; ADDRESS OF DATA STRUCTURE (OFFSET, THEN SEGMENT) 

OOOE 

37 

ZV DD ? ! ADDRESS OF ZV (OFFSET, THEN SEGMENT) 

0012 

S3 

OQ 

XV DD ? ; ADDRESS OF XV (OFFSET, THEN SEGMENT) 

— 

07 

90 

PARANS ENDS 


91 




92 



93 

? 



94 



— 

95 

DATA SEGMENT PUBLIC 'DATA' 


96 



0000 ???? 

97 

s 

TIvSAVE DW ? 


98 



— 

99 

DATA ENDS 


100 




101 

I ‘ 



102 

* 



103 



— 

104 

CODE SEGMENT PUBLIC 'CODE' 


105 
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LOC OBJ LIME SOURCE 


0000 

106 

FARCODE PRGC 

FAR 


107 




108 

10 ? 


0000 

110 

RFUN2L: 




0000 IE 

111 

PUSH 

BS 

5 SAVE REGISTERS 

0001 55 

112 

PUSH 

BP 



0002 06 

113 

PUSH 

ES 



0003 8 BEC 

114 

MOV 

BPi SP 



0005 BS 

R 115 

MOV 

AX, BATA 



0003 8EC0 

116 

MOV 

ES, AX 



OOOA C55E12 

117 

LBS 

BX, [BP3.XV 

; GET XV ADDRESS AND SEGMENT 

GOOD 9BD9Q7 

113 

FLB 

DWORD PTR EBX3 

iPUT XV ON TOP OF 8087 STACK 

0010 C55E0E 

11 ? 

LD3 

BX, CBP3.ZV 

;GET 2V ADDRESS AND SEGMENT 

0013 9BD?07 

120 

FLB 

DWORD PTR [BX] 

!PUT 2V GN TOP OF 3037 STACK 

0016 C55E0A 

121 

LBS 

BX, [BP3.DADDR 

7 GET 

START ADDRESS OF DATA STRUCTURE 

001? 8BB3 

122 

MOV 

BX, BX 

? SAVE BX 

0*018 8B0F 

123 

MOV 

CX» CBX3-XPTR 

;CX 

= XPTR 

00 IB SB6F04 

124 

MOV 

BP, CBX3.NXPT3 

5 BP 

= NXPTS 

0020 8B7F02 

125 

MOV 

Dh CBX3.ZPTR 

;DI 

= ZPTR 

0023 D1E7 

126 

r-iH 

•_>nL 

DIi 1 



0025 D1E7 

127 

SHL 

DIi 1 

?DI 

= BYTE INDEX = ELEMENT INDEX * 4 

0027 8B7706 

123 

MOV 

SI, [RX 3. NZPTS 

;si 

= NZPTS 

002A B1E6 

12 ? 

SHL 

Sit 1 



002C B1E6 

130 

SHL 

SI, 1 

?SI 

= NO. OF BYTES = NZPTS * 4 


131 

002E 132 GET.ZJKDEX: 


002E 83C308 


133 

ADD 

BX, 8 ; 

BX = START OF Z ARRAY 

0031 ?BB3i 1 


134 +2 

FCOM 

DWORD PTR [BX + DI3 ? 

COMPARE ARRAY ELEMENT WITH ZV 

0034 9B26BD3E0000 

R 

135 +1 

FSTSt-J 

ESiSTWSAVE 


003A 9B 


136 +1 

FHAIT 



003B 26A10000 

R 

137 +1 

MOV 

AX, ESiSTWSAVE 


003F 250041 


133 +1 

AND 

AX, 4100H ; 

MASK-IN COMPARISON BITS 

0042 3D0001 


13? +2 

CMP 

AX, 0100H 


0045 7434 


140 +2 

JE 

LEFT-OF-Z.START 

; EXIT IF ZV < ARRAY ELE 



HENT 






141 +2 






142 


' 

IF Z > DATA (ORIGINAL INDEX) 

0047 


143 RIGHT J3F.Z.STARTI 

*** THEN BEGIN *** 

0047 83C704 


144 

ABB 

DIi 4 

INDEX = INDEX + 4 

004A 3BFE 


145 RZSis 

CMP 

DIi SI 

WHILE (INDEX <= NO. OF BYTES) 

0040 741B 


146 

JE 

RZS2 




147 



*** DO BEGIN *** 

004E 9BD311 


143 +2 

FCOM 

DWORD PTR [BX + DI3 

COMPARE ARRAY ELEMENT WITH ZV 

0051 9E26DD3EOOOO 

R 

14? +1 

FSTSW 

ESiSTWSAVE 


0057 9B 


150 +1 

FWAIT 



0053 26A10000 

R 

151 +1 

MOV 

AX, ESiSTWSAVE 


005C 250041 


152 +1 

AND 

AX, 4100H 

; MASK-IN COMPARISON BITS 

005F 3B0001 


153 +2 

CMP 

AX, 0100H 


0062 7405 


154 +2 

JE 

RZS2 

I EXIT IF ZV < ARRAY ELEMENT 



155 +2 




0064 S3C704 


156 

ABB 

DIi 4 

; INDEX = INDEX + 4 

0067 EBE1 


157 

JMP 

RZS1 

? *** END (WHILE) *** 

0069 8 BC 6 


153 RZS2: 

MOV 

AX, SI 

5 

006 B 2B0400 


15? 

SUB 

AX, 4 

» 
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LC€ OBJ 


LIME 

SOURCE 



006E 83EF04 


160 

SUB 

DI, 4 


0071 3EF8 


161 

CMP 

DIi AX 

IF INDEX O MO. OF BYTES 

0073 7542 


162 

JNE 

FOUND-LOW. Z - 1 NBE X 

THEN FOUND_LOtLINDEX 

0075 9BDDD8 


163 

FSTP 

ST (0> 

ELSE OUT-OF-RANGE 

0078 EB4B90 


164 

JMP 

GET- X_ INDEX 




165 



*** END (IF) m 



166 



IF Z < DATA (ORIGINAL INDEX) 

007B 


167 

LEFT-0F-Z-3TART 

: 

m* THEN BEGIN m 

007B 83EF04 


168 

SUB 

DI, 4 

INDEX = INDEX - 4 

007E 83FF00 


169 

LZSl: CMP 

DI, 0 

WHILE (INDEX > 0) 

008! 7C20 


170 

JL 

LZS2 




171 



£** DO BEGIN *** 

0083 9BD311 


172 +2 

FCOM 

DWORD PTR CBX + DI) 

COMPARE ARRAY ELEMENT WITH ZV 

0086 9B26BD3E0000 

R 

173 +1 

FSTSW 

ESsSTVfSftVE 


OOSC ?B 


174 +1 

FNAIT 



0030 26A10000 

R 

175 +1 

MOV 

AX, ES: STHSAVE 


0091 250041 


176 +1 

AND 

AX, 4100H 

MASK-IN COMPARISON BITS 

0094 3B0000 


177 +2 

CMP 

AX, OOOOH 

EXIT IF ZV > ARRAY ELEMENT 

0097 740A 


178 +2 

JE 

LZS2 


0099 3B0040 


179 -f-2 

CMP 

AX, 4000H 


009C 7405 


ISO +2 

JE 

LZS2 

EXIT IF ZV = ARRAY ELEMENT 

009E 32EF04 


181 

SUB 

DI, 4 

INDEX = INDEX - 4 

00 A 1 EBDB 


182 

JMP 

LZSl 

*** END (WHILE) *** 

00A3 S3FF00 


183 

LZS2: CMP 

DI, 0 

IF INDEX O 0 

00A6 7D0F 


134 

JGE 

FOUND-LOH-Z-INDEX 

THEN FOUND-LOW-INDEX 

00A3 83C704 


185 

ADD 

DI, 4 

ELSE OUT-Cf -RANGE 

OOAB 9BBBD3 


186 

FSTP 

ST (0) ! 


OOAE 9BD9EE 


137 

FLDZ 


REPLACE TOP WITH 0 

GOBI 9BD9C9 


188 

FXCH 

ST { 1 > ' 

PUT XV ON TOP OF 8087 ST. 

00B4 EB1190 


1S9 

JMP 

GET-X-INDEX ? 




190 


i 

*** END (IF) *** 



191 




OOB7 


192 

FOUNB-LOW_Z_I NDEX : 


00B7 9BD901 


193 

FLD 

DWORD PTR EBX+DI3 

1 STORE ZL ON STACK TOP 

OOBA 9BDCE9 


194 

FSUB 

ST ( 1 ) , ST 

i (ZV-ZL) 2ND FROM TOP 

OOBD 9BD36904 


195 

FSUBR 

DWORD PTR EBX+BI+4] 

KZH-ZL) ON STACK TOP 

00C1 9BDEF9 


196 

FDIV 


l ZV-ZL) /(ZH-ZL) LEFT ON STACK 

00C4 9PB9C9 


197 

FXCH 

ST ( 1 > ;PUT XV ON TOP OF 8087 STACK 



198 




0QC7 


199 

GET_ X_ INDEX: 



00C7 03 DE 


200 

ADD 

BX, SI ;BX = START OF X ARRAY 

00C9 87F9 


201 

XCHG 

DI, CX ;DI = XPTR : CX = LOW.Z.INDEX 

OOCB B1E7 


202 

SHL 

DI, 1 


OOCB D1E7 


203 

SHL 

DI, 1 

DI = BYTE INDEX = XPTR * 4 

OOCF S7F5 


204 

XCHG 

SI, BP 

SI = NXPTS 

00B1 DIES 


205 

SHL 

SI, 1 


OOD3 B1E6 


206 

SHL 

SI, 1 

SI = NO. OF X BYTES = NXPTS * 4 

00B5 9BD811 


207 +2 

FCOM 

DWORD PTR [BX + DI] 

COMPARE ARRAY ELEMENT KITH XV 

00D8 9B26DB3E0000 

R 

208 +1 

FSTSW 

ES: STHSAVE 


CODE ?B 


209 +1 

FWAIT 



OODF 26A10000 

R 

210 +1 

MOV 

AX, ES5STWSAVE 


00E3 250041 


211 +1 

AMD 

ax, mm 

MASK-IN COMPARISON BITS 

00E6 3B0001 


212 +2 

CMP 

AX, 0100H 


00E9 7434 


213 +2 

JE 

LEFT-OF-X-START 

5 EXIT IF XV < ARRAY ELE 


KENT 
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LOC OBJ 


LINE SOURCE 


OOEB 


214 +2 

215 

216 

RIGHT-OF.X-START: 

OOEB 83C704 


217 

ADD 01, 4 

OGEE 3BFE 


218 

RXSls CMP DIi SI 

OC'FO 741B 


219 

JE RXS2 

OOF2 9BDS11 


220 

221 +2 

FCOM DWORD PTR EBX + DI3 

CK>F5 9326BB3EOOOO 

R 

222 +1 

FSTSW ESiSTWSAVE 

OOFB SB 


223 +1 

FWAIT 

OOFC 26 A 10000 

R 

224 +1 

MOV AX, ESiSTWSAVE 

0100 250041 


225 +1 

AND AX, 4100H 

0103 3D0001 


226 +2 

CMP AX, 0100H 

0106 7405 


227 +2 

JE RXS2 

0108 S3C704 


228 r2 
22? 

ADD 01, 4 

010B EBE1 


230 

JMP RXS1 

010D 8BC6 


231 

RXS2: MOV AX, SI 

01 OF 2B0400 


232 

SUB AX, 4 

0112 83EF04 


233 

SUB 01, 4 

0115 3BF3 


234 

CMP 01, AX 

0117 753F 


235 

JNE FOUND-LGW-X- INDEX 

Oil? 9BDDD8 


236 

FSTP ST(0) 

011C EB4790 


237 

JMP GET-YL 

01 IF 


ooq 

23? 

240 

LEFT-OF-X-START: 

01 IF S3EF04 


241 

SUB 01, 4 

0122 S3FF00 


242 

LXSl: CMP 01, 0 

0125 7C20 


243 

JL LXS2 

0127 9BD311 


244 

245 +2 

FCOM DWORD PTR (BX + 013 

012A 9B26DD2E0000 

R 

246 +1 

FSTSW ESiSTWSAVE 

0130 ?B 


247 +1 

FWAIT 

0131 26A10000 

R 

243 +1 

MOV AX, ESiSTWSAVE 

0135 250041 


24? +1 

AMD AX, 4100H 

0133 3D0000 


250 +2 

CMP AX, OOOOH 

013B 740A 


251 +2 

JE LXS2 

01 3D 3D0040 


252 +2 

CM? AX, 4000H 

0140 7405 


253 +2 

JE LXS2 

0142 33EF04 


254 

SUB DI, 4 

0145 EEDB 


255 

JMP LXSl 

0147 32FF00 


256 

LXS2S CMP DI, 0 

01 4 A 7D0C 


257 

JGE FOUND-LOHLX-INDEX 

CMC 830704 


253 

ADD DI, 4 

01 4F 9BDDD3 


25? 

FSTP ST(O) 

0152 9BD9EE 


260 

FLD2 

0155 EB0E90 


261 

JMP GET-YL 

0153 


262 

263 

264 

FOUND-LOW_X_INDEXi 

0158 9BD901 


265 

FLD DWORD PTR (BX+DI3 

015B 9BDCE9 


266 

FSU3 ST ( 1 ) , ST 

015E 9BD36904 


267 

FSUBR DWORD PTR [BX+DI+43 

0162 9BDEF9 


263 

FDIV 


IF X > DATAtORIGINAL INDEX) 

*** THEN BEGIN *** 

INDEX = INDEX + 4 

WHILE (INDEX O NO. OF BYTES) 

m * DO BEGIN *+* 

COMPARE ARRAY ELEMENT WITH XV 


MASK-IN COMPARISON BITS 

EXIT IF XV < ARRAY ELEMENT 

INDEX = INDEX + 4 
*** END (WHILE) *** 


IF INDEX O NO. Cf BYTES 
THEN FOUND-LOW. INDEX 
ELSE OUT_OF-RANGE 

*** END (IF) «* 

IF X < DATACORIGINAL INDEX) 

*** THEN BEGIN *** 

INDEX = INDEX - 4 
WHILE (INDEX > 0) 

DO BEGIN *** 

COMPARE ARRAY ELEMENT WITH XV 


MASK-IN COMPARISON BITS 
EXIT IF XV > ARRAY ELEMENT 


EXIT IF XV = ARRAY ELEMENT 
INDEX = INDEX - 4 
*** END (WHILE) *** 

IF INDEX O 0 

THEN FOUND-LOW-INDEX 
ELSE 0UT_0F-RANGE 

REPLACE TOP WITH 0 

*** END (IF) *** 


5 STORE XL ON STACK TOP 
5CXV-XL) ON STACK BOTTOM 
»(XH-XL) ON STACK TOP 
;(XV-XL)/(XH-XL) LEFT ON STACK 
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LOC OBJ 

LINE 

SOURCE 




269 




0165 

270 

GET.YL: 



0165 03DE 

271 

ADD 

BX, SI JBX = START OF YARRAY 

0167 8BC6 

272 

NOV 

AX, SI 

AX = NO. GF X BYTES 

0169 D1E8 

273 

SHR 

AX, 1 


016B D1E3 

274 

SHR 

AX, 1 

AX = NXPTS = NO. OF X BYTES / 4 

0160 F6E1 

275 

MUL 

CL 

AX = YA ROW OFFSET 


276 



THE AVGVE INSTR. LIMITS LENGTH CF ZARRAY 


277 



TO 255 ELEMENTS' 

01 6F 03D8 

278 

ADD 

ex, ax 

BX = YA ROW PTR = NXPTS * ZPTR 

0171 03DF 

279 

ADD 

BX, DI 

BX = YA POINTER 

0173 9BD907 

230 

FLD 

DWORD PTR [BX] 

STORE YA ON ST TP 

0176 SBC? 

281 

MOV 

AX, DI 


0173 050400 

232 

ADD 

AX, 4 


017E 3BC6 

283 

CMP 

AX, SI 


017D 7410 

234 

JE 

GET-YH J EXIT IF XV PAST 30UNDRY 

01 7F 9BB9C9 

285 

FXC.4 

ST(1) 

[ X 1 ?— XL ) / ( XH— XL > ON TOP, THEN YA 

0182 9BB94704 

286 

FLD 

DWORD PTR EEX + 43 

STORE YB ON ST. TOP 

0136 9BB3E2 

287 

FSUB 

ST, ST (2) 

(YB-YA) ON ST. TOP 

0189 93D8C9 

233 

FMUL 

ST, ST ( 1 ) 

( ( XV-XL ) / ( XH— XL ) ) ( Y3-YA ) ON ST TP 

0180 9BDEC2 

239 

FADDP 

ST (2), ST 5 YL 2ND FROM TOP OF STACK 


290 




01 SF 

291 

GET.YH: 



018F 8BC1 

292 

MOV 

fix, CX ; IF {Z_PO INTER TOO BIG) 

0191 050400 

293 

ADD 

AX, 4 ;*** DO BEGIN *** 

0194 3BC5 

294 

CMP 

AX, BP 


0196 750F 

295 

JNE 

GYHl 


0198 SBC7 

296 

MOV 

AX, DI 

IF ( X.PO INTER TOO BIG) 

019A 050400 

297 

ADD 

AX, 4 


019D 3BC6 

293 

CMP 

AX, SI 


019F 742A 

OOQ 

JE 

FOUND.YV 

THEN EXIT 

01A1 9BBBB3 

300 

FSTP 

ST 

ELSE DISCARD JUNK 

01A4 EB2590 

301 

JMP 

FOUND.YV 

AND THEN EXIT; 


302 



END (IF3 ### 

01A7 03BE 

303 

GYHl: ADD 

BX, SI 

BX = YC POINTER 

01 A9 9BD907 

304 

FLD 

DWORD PTR [BX] 

STORE YC ON STACK TOP 

02 AC 8BC7 

305 

MOV 

AX, DI 


01 AE 050400 

306 

ADD 

AX, 4 


01B1 3BC6 

307 

CMP 

AX, SI 


01 B3 740B 

308 

JE 

GET.YV ;EXIT IF XV PAST EOUNDRY 

01B5 9BB94704 

309 

FLD 

DWORD PTR [BX + 43 

'STORE YD ON STACK TOP 

01B9 9BD3E1 

310 

FSUB 

ST, ST ( 1 ) 

YD-YC ON ST TOP 

01 BC 9BDECA 

311 

FKULP 

ST (2) * ST 

( ( XV-XL ) / ( XH-XL ) ) ( YD-YC ) 2ND FROM TOP 

01 BF 9BBEC1 

312 

FADD 

5YH ON TOP OF 8087 


313 




01C2 

314 

GET.YV: 

?YH Tip, THEN YL, THEN SLOPE 

0102 9BB3E1 

315 

FSUB 

ST, ST ( 1 > ’YH-YL ON ST. TOP 

01 C5 9BDECA 

316 

FMULP 

ST(2)» ST 

( ( Z V-ZL ) / ( ZH-ZL ) ) ( YH-YL ) 2ND FROM TOP 

01C8 9BBEC1 

317 

FADD 

;yv ON ST. TOP 


318 




01CB 

319 

FOUND. YV i 



01CB 

320 

SAVE.NEW.PO INTERS: 


01 CB 3BBA 

321 

MOV 

BX, DX 


01 CD D1EF 

322 

SHR 

DI, 1 


01 CF D1EF 

323 

SHR 

01, 1 



31 
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LOG OBJ LINE SOURCE 


01D1 S73F 

324 

MOV 

CBX3.XPTR, DI 

, SAVE XPTR 

01D3 DIE? 

325 

SHR 

CX, 1 


01D5 DIE? 

326 

SKR 

CX, 1 


01D7 8?4F02 

327 

ooo 

MOV 

[BX3.ZPTR, CX 

; SAVE ZPTR 

01DA 

wi_L‘ 

32? 

EPILOGUE: 



OiBA 07 

330 

POP 

ES 


01DB 5D 

331 

POP 

BP 

•RESTORE REGISTER' 

01 DC IF 

332 

POP 

DS 


01DD CAOCOO 

333 

RET 

12 



334 

335 FARCODE ENDP 

336 

337 CODE ENDS 

333 

340 •****m***#*mm*MM#m#***M^m*mM*M*m«*#fr**#*#m^iHH*’ 

341 

342 END 
ASSEMBLY COMPLETE* NO ERRORS FOUND 
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SERIES-III 8086/37/88/186 MACRO ASSEMBLER V2.0 ASSEMBLY OF MODULE RFUH3S 
OBJECT MODULE PLACED I N :F1:RFUN3S.0BJ 
ASSEMBLER INVOKED BY: ASM 86.86 :F1:RFUN3S.SRC 


LOC OBJ 


LHC SOURCE 


1 

2 

3 

4 

5 

6 
7 
3 

9 




* ROUTINE NAME: RFUN3S * 

* DATE: SEP. 6 , 1933 * 

* PROGRAMMER'S NAME: MICHAEL MACKIN * 

* PURPOSE: TO INTERPOLATE A SIMPLE * 

* FUNCTION FOR SMALL-MODEL * 

* PROGRAMS 4 




10 

11 

12 

13 

14 

15 

16 
17 
13 


THIS SUBROUTINE OUTPUTS THE ESTIMATED VALLE OF A FUNCTION YV = F(XV?ZV) 

BY INTERPOLATING FROM KNOW SOLUTION POINTS? USING THE RELATIONS 

XF = ( ( ZV-ZL > / ( ZH-ZL ) H XB-XC ) + XC 

XG = ( ( ZV-ZL > /( ZH-ZL > M XE-XD I + XD 

YF = ( (ZV-ZL)/ ( ZH-ZL ) ) (YB-YC) + YC 

YG = ( ( ZV-ZL ) / ( ZH-ZL ) ) ( YE-YD > + YD 

YV = ( ( XV-XF ) / ( XG-XF ) ) ( YG-YF ) + YF 


INPUT REQUIREMENTS: 


19 

20 

21 


1) THE FOLLOWING ADDRESSES HILL BE PASSED TO THIS SUBROUTINE: 

A. ON TOP OF THE 8086 STACK: 

THE ADDRESS (OFFSET) OF A DATA AREA WITH THE FOLLOWING FORMAT 




22 

23 

BSTRUC 


24 


0000 

25 

XPTR 

0002 

26 

ZPTR 

0004 

27 

NXPTS 

0006 

23 

NZPTS 


29 ; 

ZARRAY 


30 ! 

XARRAY 


31 s 

YARRAY 


32 


— 

33 

DSTRUC 


STRUG 



ON 

9 

5X ARRAY INDEX 

DW 

9 

JZ ARRAY INDEX 

DW 

7 

5 NO. OF X BREAKPOINTS ON EACH Z LINE 

DW 

7 

SNO. OF ELEMENTS IN Z ARRAY 

DD 

NZPTS DUP(?) 

?Z ARRAY (255 ELEMENTS MAXIMUM) 

DO 

NXPTS DUPf?> 

?X ARRAY 

DD 

NYPTS DUP(?) 

?Y ARRAY 

ENDS 




34 


35 

36 


B. ON THE 8037 STACK: ON TOP ZV, THEN XV 


37 

38 
3? 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


2) IT IS REQUIRED THAT THE 8087 CHIP HAVE 5 EMPTY REGISTERS 
WHEN THIS INTERPOLATION ROUTINE IS CALLED. 

3) NOTE: NO TWO ADJACENT ELEMENTS OF X ARRAY? YARRAY, OR ZARRAY 
MAY HAVE IDENTICAL ELEMENTS. THIS WILL RESULT IN DIVISION 
BY ZERO. 

OUTPUT EFFECTS: 

1) REGISTERS DESTROYED: AX? BX? CX? DX? DI? SI 

2) INTERPOLATED RESULT RETURNED ON TOP OF 8087 STACK 

3) XPTR LOCATION UPDATED TO INDEX X VALLE JUST PRECEDING XV. 
ZPTR LOCATION UPDATED TO INDEX Z VALLE JUST PRECEDING ZV. 
(TO SPEED FUTURE SEARCHES FOR X AND Z) 
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LOC OBJ LINE SOURCE 



51 ? 


EXAMPLE PL/M CALL 


52 5 




53 3 


PLHPGMs DO; 


54 




55 5 


DECLARE DSTRUC STRUCTURE ( 


56 


X INDEX INTEGER, Z INDEX INTEGER, NXPTS INTEGER, NZPTS INTEGER 


57 


ZARRAY<3) REAL, X ARRAY 05) REAL, YARRAY(15) REAL) PUBLIC 


58 


INITIAL ( 


59 


/* XINDEX */ 0, 


60 


/* ZINDEX */ 0, 


61 


/# NXPTS */ 5, 


62 


/* NZPTS */ 3, 


63 


/* ZARRAY */ 0.0, 10.0, 20.0, 


64 


/* X ARRAY */ 0.0, 2.0, 4.0, 6.0, 3.0, 


65 


1.0, 3.0, 5.0, 7.9, 9.0, 


66 


1.5, 1.75, 1.3, 5.0, 10.0, 


67 


/* YARRAY */ 0.0, 2.0, 4.0, 6.0, 8.0, 


68 


10.0, 12.0, 14.0, 16.0, 18.0, 


69 


20.0, 22.0, 24.0, 26.0, 28.0 )5 


70 




71 


RFUN3SS PROCEDURE ( X , Z , STRUC.ADDR ) REAL EXTERNAL; 


72 


DECLARE X REAL; 


73 


DECLARE Z REAL; 


74 


DECLARE STRUC.ADDR POINTER? 


75 


END RFUN3S; 


76 




77 


R3TEST; PROCEDURE PUBLIC? 


73 


DECLARE (X,Z,Y) REAL; 


79 


X = 3.2; 


90 


z = 12.0; 


31 


Y = RFUN3S ( X , Z i 3DSTRUC ) ? 


82 


CALL PRINTOUT (Y)? 


33 


END; 


84 


END R3TEST? 


35 




86 


END PLMPGM; 


87 




S3 



89 




90 




91 


NAME RFUN3S 


92 

CGRGUP 

GRC4JP CODE 


93 

DGROUP 

GROUP DATA 


94 


ASSUME CS! CGRGUP, DS? DGROUP 


95 


PUBLIC RFUN3S 


96 




97 

;^****#«****#«##Mmwm***M***M*M**MMmmM***s**s*m*** 


98 




99 

5 

STACK STRUCTURE 


100 



— 

101 

PARAMS 

STRUC 


102 



0000 

103 

OLB-BP 

DU ? ? SAVED BP REGISTER 

0002 

104 

RETURN 

DU ? 5 RETURN ADDRESS 

0004 

105 

BABBR 

DU ? ; ADDRESS OF DATA STRUCTURE 
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LOC OBJ 


LINE 

SOURCE 





106 




— 


107 

PAR AMS ENDS 





103 






10? 




110 




111 




— 


112 

BATA SEGMENT PUBLIC 7 DATA' 




113 




0000 ???? 


114 

STHSAVE BN 

7 




115 




— 


116 

DATA ENDS 





117 






118 

; ************** 

*************************** 




ii? 

120 

'wwmwwzwwww************************************************** 






— 


121 

CODE SEGMENT PUBLIC 'CODE' 




122 




0000 


123 

124 

125 

NEARCOBE 

PROC NEAR 




126 




0000 


127 

RFLIN3S: 





128 


’SAVE REGISTERS 

0000 55 


12? 

PUSH 

PP 


0001 3BEC 


130 

MOV 

BP, SP 


0003 8B5E04 


131 

MOV 

BX, CBP3.DADDR ;GET START ADE*RESS OF DATA STRUCTURE 

0006 SBD3 


132 

MOV 

DX, BX 

SAVE BX 

0003 8B2F 


133 

MOV 

BP, EBX3.XPTR 

BP = XPTR 

OOOA 8B4F04 


134 

MOV 

CX, EBX3.NXPTS 

CX = NXPTS 

OOOB 8B7F02 


135 

MOV 

DIi EBX3.ZPTR ;DI = ZPTR 

0010 B1E7 


136 

SHL 

DI, 1 


0012 B1E7 


137 

SHL 

DIi 1 

DI = BYTE INDEX = ELEMENT INDEX * 4 

0014 8B7706 


133 

MOV 

SI, CBX3.NZPTS 5 SI = NZPTS 

0017 D1E6 


13? 

SHL 

SI, 1 


0019 B1E6 


140 

SHL 

SI, 1 5 SI = NO. OF BYTES = NZPTS * 4 



141 




00 IB 


142 

GET.ZJNBEX: 



001 B 83C303 


143 

ADD 

BX, 8 ;BX = START OF Z ARRAY 

00 IE 9BD811 


144 +2 

FCOM 

DWORD PTR EBX + DI3 

COMPARE ARRAY ELEMENT WITH ZV 

0021 9BDD3E0000 

R 

145 +1 

FSTSW 

STHSAVE 


0026 9B 


146 +1 

FKAIT 



0027 A10000 

R 

147 +1 

MOV 

AX, STHSAVE 


002A 250041 


143 +1 

AND 

AX, 4100H 

MASK-IN COMPARISON BITS 

002D 3D0001 


14? +2 

CMP 

AX, OiOOH 


0030 7432 


150 +2 

JE 

LEFT-OF_Z_START 

? EXIT IF ZV < ARRAY 




HENT 





151 +2 

152 



IF Z > DATAtORIGIKflL INDEX) 

0032 


153 

RIGHT J3F.Z.START: 

*** THEN BEGIN *** 

0032 83C704 


154 

ADD 

DI, 4 

INDEX = INDEX + 4 

0035 3BFE 


155 

RZSl: CMP 

DI, SI 

WHILE (INDEX <= NO. OF BYTES) 

0037 741? 


156 

JE 

RZS2 




157 



*** DO BEGIN 

003? 9BD811 


153 +2 

FCOM 

DWORD PTR EBX + DI3 

COMPARE ARRAY ELEMENT WITH ZV 

003 C 9BBB3E0000 

R 

15? +1 

FSTSW 

STHSAVE 
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LOC OBJ LINE SOURCE 


0041 9B 


160 +1 


FWAIT 



0042 A 10000 

R 

161 +1 


MOV 

AX, STHSAVE 


0045 250041 


162 +1 


AND 

AX, 4100H ; 

MASK-IN COMPARISON BITS 

0043 3D0001 


163 r2 


CMP 

AX, 0100H 


004B 7405 


164 +2 


JE 

RZS2 

EXIT IF ZV < ARRAY ELEMENT 



165 +2 





004D 33C704 


166 


ADD 

DI. 4 

INDEX = INDEX + 4 

0050 EBE3 


167 


JHP 

RZSi 

**# END (WHILE) *** 

0052 8BC6 


163 

RZS2: 

MOV 

AX, SI 


0054 2D0400 


169 


SUB 

AX, 4 


0057 83EF04 


170 


SUB 

DI, 4 


005A 3EF3 


171 


CNF 

DI, AX 

IF INDEX O NO. OF BYTES 

005C 7540 


172 


JNE 

FOUND -LOW_Z_ I NDEX 

THEN FOUND-LOW-INDEX 

005E 9BDDD3 


173 


FSTP 

ST(O) 

ELSE OUT-OF-RANGE 

0061 EB4B90 


174 


JHP 

GET-X- INDEX 




175 




END (IF) 



176 




IF Z < DATA(ORIGINAL INDEX) 

0064 


177 

LEFT_OF_Z-START: 

*** THEN BEGIN *** 

0064 83EF04 


173 


SUB 

DI, 4 

INDEX = INDEX - 4 

0067 S3FF00 


179 

LZSl: 

CMP 

DI, 0 

WHILE (INDEX > 0) 

006A 7C1E 


180 


JL 

LZS2 




181 




m * DO BEGIN *** 

006C 9BD311 


182 +2 


FCOM 

DWORD PTR EBX + DI] 

COMPARE ARRAY ELEMENT WITH ZV 

006F 9BDD3E0000 

R 

183 +1 


FSTSW 

STHSAVE 


0074 9B 


184 +1 


FWAIT 



0075 A10000 

R 

135 +1 


MOV 

AX, STHSAVE 


0073 250041 


186 +1 


AND 

AX, 41 COM 

MASK-IN COMPARISON BITS 

007B 3D0000 


137 +2 


CMP 

AX, OOOOH 

EXIT IF ZV > ARRAY ELEMENT 

007E 740A 


183 +2 


JE 

LZS2 


0080 3D0040 


18? +2 


CMP 

AX, 4000H 


0033 7405 


190 +2 


JE 

LZS2 

EXIT IF ZV = ARRAY ELEMENT 

0035 83EF04 


191 


SUB 

DI, 4 

INDEX = INDEX - 4 

0088 EBDD 


192 


JMP 

LZSl 

*** END (WHILE) «* 

OOSA 33FF00 


193 

LZ32: 

CMP 

DI, 0 

IF INDEX O 0 

008D 7B0F 


194 


JOE 

FOUND-LOO - 1 NDEX 

THEN FOUND-LOW-INDEX 

OOSF 83C704 


195 


ADD 

DI, 4 

ELSE OUTJF-RANGE 

0092 9BBBD3 


196 


FSTP 

ST (0) 


0095 9BD9EE 


197 


FLDZ 


REPLACE TOP WITH 0? 

0098 9BB9C9 


193 


FXCH 

STU) 


009 B EB1190 


19? 


JMP 

GETJLINBEX 




200 




END (IF) 

009E 


iOl 

202 

FOUND-LOW. Z. I NDEX : 


009E 9BD901 


203 


FLD 

DWORD PTR [BX+DI3 

•STORE ZL ON STACK TOP 

00 A 1 9BBCE9 


204 


FSUB 

ST ( 1 ) , ST 

5 (ZV-ZL) 2ND FROM TOP 

00A4 9B086904 


205 


FSUBR 

DWORD PTR [BX+DI+43 

>(ZH-ZL) ON STACK TOP 

WAS 9BBEF9 


206 


FDIV 


;( ZV-ZL )/(ZH-ZL) LEFT ON STACK 

OOAB 9BD9C9 


207 


FXCH 

ST ( 1 > 

?PUT XV ON TOP OF 8087 STACK 



203 





00 AE 


209 

OET-X. 

.INDEX: 



OOAE 03DE 


210 


ADD 

BX, SI 

?BX = START OF XARRAY 

OOBO 3BC7 


211 


MOV 

AX, DI 

,*AX = NO. OF Z BYTES 

00B2 F6E1 


212 


MIL 

a 

?AX = XC ROW OFFSET=N. Cf ZBYTES*NXPTS 

00B4 03D3 


213 


ADD 

BX, AX 

;BX = XC ROW PTR 

00B6 37FB 


214 


XCHG 

DI, BP 

5DI = XPTR : BP = LOH-Z-INDEX 


36 



8036/87/88/186 MACRO ASSEMBLER RFUN3S 


08/30/84 PAGE 5 


LOG OBJ 


LIKE SOURCE 


OOB3 D1E7 
OOBA B1E7 
OOBC 87F1 
OOBE 01E6 
OOCO B1E6 
00C2 9BD901 
00C5 8BC5 
0007 050400 
OOCA 3BC1 
OOCC 7410 
OOCE 03BF 
OOBO 9BD900 
00B3 2BDF 
00B5 9BDSE1 
00B3 9BB8CB 
OODB 9BBEC1 
CODE 9BD8D1 
GOBI 9BDD3E0000 
00E6 9B 
00E7 A 10000 
OOEA 250041 
COED 3D0000 
OOFO 745A 


215 

216 

217 

218 

219 

220 
223 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 +2 

232 +1 

233 +1 

234 +1 

235 +1 

236 +2 

237 +2 

238 +2 


GXls 


SHL 

DI, 1 


SHL 

DI, 1 

!DI = BYTE INDEX = XPTR * 4 

XCHG 

SI. CX 

!SI = NXPTS : CX = NO. OF Z BYTES 

J 

SHL 

SI, 1 

SHL 

SI, 1 

5 SI = NO. OF X BYTES = NXPTS » 4 

FLD 

DWORD PTR IBX + DI) 

! STORE XC ON TO? OF STACK 

MOV 

AX, BP 


ADD 

AX, 4 


CMP 

AX, CX 


JE 

GX1 

5EXIT IF Z TOO HIGH 

ADD 

BX, DI 


FLD 

DWORD PTR IBX + SI) 

i STORE XB ON TOP OF STACK 

SUB 

BX, DI 


FSUB 

ST, STU) 

iXB-XC ON ST. TO? 

FflUL 

ST, ST ( 3 > 

» ( ( ZV-ZL )/( ZH-ZL) ) ( XB-XC) ON ST. TOP 

FADD 


!XF ON ST. TOP 

FCOM 

ST(!) 

? COMPARE ARRAY ELEMENT WITH XV 

FSTSW 

STHSAVE 

FHAIT 



MOV 

AX, STHSAVE 


AND 

AX, 4100H 

! MASK-IN COMPARISON BITS 

CMP 

AX, OOOOH 

; EXIT IF XV > ARRAY ELEMENT 

JE 

LEFT.OF.X.START 





239 




00F2 


240 

RIGHT_OF_X_START: 

«X>F2 83C704 


241 


ADD 

DI, 4 

OOFS 3BFE 


242 

RXS1 : 

CMP 

DI, SI 

OOF7 743B 


243 


JE 

RXS2 



244 




OOF? 9BD901 


245 


FLD 

DWORD PTR EBX + DI] 

OOFC 8BC5 


246 


MOV 

AX, BP 

OOFE 050400 


247 


ADD 

AX j 4 

010! 3BC1 


243 


CMP 

AX, CX 

0103 7410 


249 


JE 

RXS3 

0105 03DF 


250 


ADD 

BX, DI 

0107 9BD900 


251 


FLD 

DWORD PTR [BX + SI] 

010A 2BDF 


252 


SUB 

BX, DI 

010C 9BD8E1 


253 


FSUB 

ST, STU) 

010F 9BDSCC 


254 


FMUL 

ST, ST(4) 

0112 9BDEC1 


255 


FADD 


0115 9BD8D2 


256 +2 

RXS3: 

FCOM 

ST(2) 

0118 9BDD3E0000 

R 

257 +1 


FSTSW 

STHSAVE 

01 ID 93 


253 +1 


mu 


01 IE A 10000 

R 

259 +1 


MOV 

AX, STHSAVE 

0121 250041 


260 +1 


AND 

AX, 4200K 

0124 3D0000 


261 +2 


CMP 

AX, OOOOH 

0127 740B 


262 +2 


JE 

RXS2 



263 +2 




0129 9BD9C9 


264 


FXCH 

STU) 

012C 9BDBB3 


265 


FSTP 

ST(O) 

012F 83C704 


266 


ADD 

DI, 4 

0132 EBC1 


267 


JMP 

RXS1 

0134 9BD9C9 


268 

RXS2: 

FXCH 

STU) 

0137 8BC6 


269 


MOV 

AX, SI 


IF X > DATAIORIGINAL INDEX) 

*** THEN BEGIN *** 

INDEX = INDEX + 4 

HHILE (INDEX <= NO. OF BYTES) 

*** DO BEGIN **# 

STORE XD ON TOP OF STACK 


EXIT IF Z TOT HIGH 
STORE XE ON TOP OF STACK 
XE-XD ON ST. TOP 

( ( 2V-ZL ) / ( ZH-ZL) } ( XD-XE) ON ST. TOP 
XG ON ST. TOP 

COMPARE ARRAY ELEMENT WITH XV 


MASK-IN COMPARISON BITS 
EXIT IF XV > ARRAY ELEMENT 


XF = XG 

DISCARD OLD XF 
INDEX = INDEX + 4 
*** END (HHILE) *** 


37 
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LOC OBJ 

LINE 

SOURCE 




0139 2D0400 

270 


SUB 

AX? 4 


013C 83EF04 

271 


SUB 

DI, 4 


013F 3BF8 

272 


CMP 

QIi AX 

IF INDEX G NO. Cf BYTES 

OH! 7565 

273 


JNE 

FOUND-LOO- 1 NDEX 

THEN FOUND-LOW_INDEX 

0143 9BBBDS 

274 


FSTP 

ST(O) 

ELSE OUT-OF_RANGE 

0146 9BBBD8 

275 


FSTP 

ST (0) 

DISCARD XG, AND XV 

0149 EB7290 

276 


JMP 

GET-ANSWER 



277 




*** END (IF) *** 


278 




IF X < DATA (ORIGINAL INDEX) 

014C 

279 

LEFT.GF.X-STftRT 


*** THEN BEGIN *** 

014C 83EF04 

230 


SUB 

DI? 4 

INDEX = INDEX - 4 

014F 33FF00 

231 

lxsi; 

CMP 

DI? 0 

WHILE (INDEX > 0) 

0152 7C40 

232 


JL 

LXS2 



233 




DO BEGIN 

0154 9BB901 

284 


FLB 

DWORD PTR EBX + DI3 

STORE XD ON TOP OF STACK 

0157 8BC5 

285 


MOV 

AX, BP 


0159 050400 

286 


ADD 

AX, 4 


015C 3BC1 

287 


CMP 

AX, CX 


015E 7410 

283 


JE 

LXS3 

EXIT IF Z TOO HIGH 

0160 03BF 

239 


ADD 

BX, DI 


0162 9BD900 

290 


FLD 

DWORD PTR CBX + SI 3 

STORE XE ON TOP OF STACK 

0165 2BBF 

291 


SUB 

BX, DI 


0167 9BD3E1 

292 


FSUB 

ST, ST ( 1 ) 

XB-XC ON ST. TOP 

016A 9BB8CC 

293 


FMUL 

ST, ST (4) 

( ( ZV-ZL ) / ( ZK-ZL > ) ( XB-XC ) ON 

016B 9BBEC1 

294 


FADD 


XG ON ST. TOP 

0170 9BB8D2 

295 +2 

LXS3: 

FCGM 

ST <2) 

COMPARE ARRAY ELEMENT WITH 

0173 9BDD3E0000 

R 296 +1 


FSTSN 

STWSAVE 


0178 9B 

297 +1 


FWAIT 



0179 A 10000 

R 298 +1 


MOV 

AX, STWSAVE 


017C 250041 

299 +1 


AND 

AX, 4100H 

5 MASK-IN COMPARISON BITS 

017F 3D0001 

300 +2 


CMP 

AX, 0100H 


0182 7410 

301 +2 


JE 

LXS2 

5 EXIT IF XV < ARRAY ELEMENT 

0184 3D0040 

302 +2 


CMP 

AX i 4000H 


0187 74 OB 

303 +2 


JE 

LXS2 

I EXIT IF XV = ARRAY ELEMENT 

0189 9BB9C9 

304 


FXCH 

ST ( 1 ) 

; XF=XG 

013C 9BBDB8 

305 


FSTP 

ST(O) 

; DISCARD OLD XF 

018F 83EF04 

306 


SUB 

DI, 4 

? INDEX = INDEX - 4 

0192 EBBB 

307 


JMP 

LXS1 

? *** END (WHILE) *** 

0194 83FF00 

303 

LXS2: 

CMP 

DI, 0 

5 IF INDEX O 0 

0197 7B0F 

309 


JGE 

FOUNB-LON_X_ INDEX 

; THEN FOUND-LOW- I NDEX 

0199 83C704 

310 


ADD 

DI, 4 

; ELSE OUT-OF-RANGE 

019C 9BBBB3 

311 


FSTP 

ST (0) 


01 9F 9BDBDS 

312 


FSTP 

ST ( 0 ) 

5 DISCARD XC 

01A2 9BB9EE 

313 


FLD2 


; REPLACE TOP WITH 0 

01A5 EB1690 

314 


JMP 

GET-ANSWER 



315 




?*** ENDCIF) #** 


316 





01 AS 

317 

FOUND -LGi-LX- 1 NBEX : 

; XF? XG, XV; Z SLOPE ON ST. 

01AS 9BBCEA 

313 


FSUB 

ST(2), ST 

? XF, XG, (XV-XF), ZSLOPE 

01AB 9BBEE9 

319 


FSUB 


; (XG-XF), (XV-XF), ZSLOPE 

01AE 9BBEF9 

320 


FDIV 


; ( (XV-XF) /(XG-XF)), ZSLOFE 

01B1 8BC5 

321 


MOV 

AX, BP 

5 IF Z NOT TOO HIGH 

01B3 050400 

322 


ADD 

AX, 4 


01B6 3BC1 

323 


CMP 

AX, CX 


01BS 7403 

324 


JE 

GET-ANSWER 



'. TOP 


38 
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LOC OBJ 

LINE 

SOURCE 



01 BA 9BD9C9 

325 

FXCH 

ST(1) 

5 THEM (ZV-ZL)/(ZH-ZL> ON ST. TOP 


326 




01BD 

327 

GET-ANSWER: 



01BB 8BC6 

32S 

MOV 

AX, SI 

5AX = NO. OF X BYTES 

01 BF DIES 

329 

SHR 

AX, 1 


01C1 DIES 

330 

SHR 

AX, 1 

5 AX = NXPTS = NO. OF X BYTES / 4 


331 



;cx = NO OF Z BYTES 

01C3 F6E1 

332 

Mil 

CL 

IAX = YC ROW OFFSET 


333 



5 THE avgve INSTR. limits length of zarray 


334 



5 TO 255 ELEMENTS 

01 C5 03D8 

335 

ADD 

BX, AX 

?BX = YC ROW PTR = NXPTS * ZPTR 

01C7 03DF 

336 

ADD 

BX, DI 

5BX = YC POINTER 

0109 9BD907 

337 

FLD 

DWORD PTR [BX] 

? STORE YC ON ST TOP 

01CC 8BC5 

333 

MOV 

AX, BP 


DICE 050400 

339 

ADD 

AX, 4 


01D1 3BC1 

340 

CMP 

AX, CX 


01D3 7510 

341 

JNE 

Z-NOT-HIGH 



342 




01 D5 

343 

Z_T0O_HIGH: 


5 IF Z OUT OF RANGE 

01D5 SBC7 

344 

MOV 

AX, DI 

;*** Dil BEGIN *** 

01D7 050400 

345 

ADD 

AX, 4 

; IF X OUT OF RANGE 

01 DA 3BC6 

346 

CMP 

AX, SI 


01 DC 7444 

347 

JE 

FOUND— YV 

? THEN EXIT 

01 DE 

343 

X.VALID: 


? CZ TOO HIGH BUT X IS IN RANGE) 

01 DE 9BD94704 

349 

FLD 

DWORD PTR [BX + 4] 

5 STORE XD 

01E2 EB3290 

350 

JMP 

GET.YV 

;*#* END CIF) m 


351 




01E5 

352 

Z-NOT-HIGH: 


5 IF Z VALID 

01E5 SBC7 

353 

MOV 

AX, DI 

;*** DO BEGIN *** 

01 E7 050400 

354 

ADD 

AX, 4 

l IF X OUT OF RANGE 

01EA 3BC6 

355 

CMP 

AX, SI 

, 

01 EC 7423 

356 

JE 

Z-VALID 

? THEN EXIT 

01 EE 

357 

BOTH. VALID: 


? ELSE (BOTH XV AND ZV ARE IN RANGE) 

01 EE 9BD9C9 

353 

FXCH 

ST( 1 ) 

;<ZV-ZL)/(ZH-ZL) Of'] TOP, THEN YC 

01F1 9BD900 

35? 

FLD 

DWORD PTR [BX + SI 3 

? STORE YB ON ST. TOP 

01F4 9BDSE2 

360 

FSUB 

ST, ST(2) 

? ( YB-YC ) ON ST. TOP 

01F7 9BDSC? 

361 

FMUL 

ST, ST ( 1 ) 

; ( ( ZV— ZL ) / ( ZH-ZL ) ) ( YB-YC } ON ST TOP 

01FA 9BDEC2 

362 

FADDP 

ST ( 2 ) , ST 

;YF 2ND FROM TOP OF STACK 

01FD 9BD94704 

363 

FLD 

DWORD PTR [BX + 43 

5 STORE YD ON STACK TOP 

0201 9BD94004 

364 

FLD 

DWORD PTR [BX + SI + 43 

? STORE YE ON STACK TOP 

0205 9BDSE1 

365 

FSUB 

ST, ST(1> 

; YE-YD ON ST TOP 

0203 9BDECA 

366 

FMULP 

ST (2) » ST 

; ((ZV-ZL)/ (ZH-ZL) MYE-YD) 2ND FROM TOP 

020B 9BDEC1 

367 

FADD 


?YG ON TOP OF 8087 

020E EP0690 

363 

JMP 

GET-YV 

?*** END (IF) *** 


369 




0211 

370 

Z-VALID: 


? (X OUT OF RANGE BUT Z VALID) 

0211 03DE 

371 

ADD 

BX, SI 


0213 9BD907 

372 

FLD 

DWORD PTR [BX3 



373 




0216 

374 

GET-YV: 


!YG TOP, THEN YF, THEN SLOPE 

0216 9BD3E1 

375 

FSUB 

ST, ST( 1 ) 

5YG-YF ON ST. TOP 

0219 9BDECA 

376 

FMULP 

ST(2), ST 

5 ( (XV-XF)/(XG-XF) ) (YG-YF) 2ND FROM TOP 

021 C 9BBEC1 

377 

FADD 


5YV ON ST. TOP 

021 F EB0190 

378 

JMP 

FGUND-YV 



37? 





39 
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LOG OBJ 

LINE 

SOURCE 




3S0 




0222 

m 

F0UMD.YV: 



0222 

332 

SAVE_NEM_POINTERS: 


0222 8BBA 

383 

MOV 

BX» DX 

J RESTORE BX TO START OF DATA STRIJC. 

0224 D1EF 

384 

SHR 

DIi 1 


0226 DIET 

335 

SHR 

DIi 1 


0228 893F 

336 

MOV 

EBXLXPTRt DI 

;SAVE XPTR 

022A DIED 

CO 

CO 

SHR 

BP, 1 


022C DIED 

383 

SHR 

BP, 1 


022E 896F02 

399 

MOV 

£BX].2PTRi BP 

; SAVE ZPTR 


390 




0231 

391 

EPILOGUE: 



0231 5D 

392 

POP 

BP 

! RESTORE REGISTERS 

0232 C20100 

393 

RET 

1 



394 





395 

NEARCODE 

EWDP 



396 




— 

397 

CODE ENDS 




398 





399 



400 



401 





402 

END 




ASSEMBLY COMPLETE, NO ERRORS FOUND 


40 
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SERIES— I I I 8086/87/88/186 MACRO ASSEMBLER V2.0 ASSEMBLY OF MODULE RFUN3L 
OBJECT MOBILE PLACED IN :Fl:RFUN3L.0BJ 
ASSEMBLER INVOKED BY: ASM86.86 :F1:RFUN3L.SRC 


LOG OBJ 


0000 

0002 

0004 

0006 


LINE SOURCE 


1 

2 

3 

4 

5 

6 

7 

8 
Q 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
2? 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
43 

49 

50 


* ROUTINE NAME: RFUM3L * 

* DATE: SEP. 8, 1933 * 

* PROGRAMMER'S NAME: MICHAEL NACKIN * 

* PURPOSE: TO INTERPOLATE A SIMPLE * 

* FUNCTION FOR LARGE-MODEL * 

* PROGRAMS * 

THIS SUBROUTINE OUTPUTS THE ESTIMATED VALUE OF A FUNCTION YV = F<XV,ZV> 

BY INTERPOLATING FROM KNOWN SOLUTION POINTS, USING THE RELATIONS 

XF = { ( ZV-ZL ) / ( ZH-ZL ) ) ( XB-XC ) + XC 

XG = ( (ZV-ZL) /( ZH-ZL ))(XE-XD) + XD 

YF = { ( ZV-ZL) / ( ZH-ZL ) ) ( YB-YC) + YC 

YG = ( ( ZV-ZL ) / ( ZH-ZL ) > ( YE-YB > + YD 

YV = ( ( XV— XF ) / ( XG-XF ) ) ( YG-YF ) + YF 

INPUT REQUIREMENTS: 

1) THE FOLLOWING ADDRESSES HILL BE PASSED TO THIS SUBROUTINE: 

A. ON TOP OF THE 3086 STACK: 

THE ADDRESS (OFFSET) OF A DATA AREA WITH THE FOLLOWING 
FORMAT 


DSTRUC 

STRUC 



XPTR 

DW 

7 

?X ARRAY INDEX 

ZPTR 

m 

7 

!Z ARRAY INDEX 

NXPTS 

DW 

7 

5N0. OF X BREAKPOINTS ON EACH Z LINE 

NZFTS 

DW 

7 

5 NO. OF ELEMENTS IN Z ARRAY 

ZARRAY 

DD 

NZPTS DUP(?) 

;Z ARRAY (255 ELEMENTS MAXIMUM) 

XARRAY 

DD 

NXPTS DUP(?) 

5X ARRAY 

YARRAY 

DD 

NYPTS DUP(?) 

?Y ARRAY 

DSTRUC 

ENDS 




B. ON THE 8037 STACK: ON TOP ZV, THEN XV 

2) IT IS REQUIRED THAT THE 8037 CHIP HAVE 5 EMPTY REGISTERS 
WHEN THIS INTERPOLATION ROUTINE IS CALLED. 

3) NOTE: NO TUG ADJACENT ELEMENTS OF X ARRAY, YARRAY, OR ZARRAY 

MAY HAVE IDENTICAL ELEMENTS. THIS WILL RESULT IN DIVISION BY ZERO. 

OUTPUT EFFECTS: 

1) REGISTERS DESTROYED: AX, BX, CX, DX, DI, SI 

2) INTERPOLATED RESULT RETURNED ON TOP OF 8087 STACK 

3) XPTR LOCATION UPDATED TO INDEX X VALUE JUST PRECEDING XV. 

ZPTR LOCATION UPDATED TO INDEX Z VALUE JUST PRECEDING ZV. 

(TO SPEED FUTURE SEARCHES FOR X AND Z) 


41 
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LQC OBJ LINE SOURCE 


0000 

0002 

0004 

0006 

OOOA 

OOOE 

0012 


0000 ???? 


51 

52 

53 

54 

55 

56 

57 
53 

59 

60 
61 
62 

63 

64 

65 

66 
67 
63 
6? 

70 

71 

72 

73 

74 

75 

76 

77 

78 
7? 
80 
81 
32 

83 

84 
35 
86 
87 
83 
8 *? 

90 

91 

92 

93 

94 

95 

96 

97 
93 
99 
100 
101 
102 

103 

104 

105 


EXAMPLE FORTRAN CALL 

SUBROUTINE F3CALL 

DIMENSION XARRAYU5), ZARRfiY<3>. YARRAYU5) 

INTEGER X INDEX i ZINDEX, NXPTS, NZPTS 

COMMON /STUFF/XINBEX,ZINDEX,NXPT$,NZPTS,ZARRAY, X ARRAY, YARRAY 
DATA XINBEXi ZINDEX, NXPTS, NZPTS/O, 0,5,3/ 

DATA ZARRAY/O. 0,10.0,20.0/ 

DATA XARRAY/Q. 0,2. 0,4. 0,6.0, 8. 0, 1., 3., 5., 7. 9,?., 1.5, 1.75, 1.8,5. ,10./ 
DATA YARRAY/0.0,2. ,4., 6. ,3., 10., 12., 14., 16., 18., 20. ,22., 24., 26., 23./ 
X=9.1 
Z=11.3 

Y=RFUN3L ( X , Z , X INDEX ) 

FRINT *, Y 
END 


NAME RFUM3L 
CGROUP GROUP CODE 

DGROUP GROUP DATA 

ASSUME CStfGRGUP, DS: DGROUP, ESs DGROUP 
PUBLIC RFUN3L 

1 STACK STRUCTURE 

PARAMS STRUC 

OLD.ES DU ? S SAVED ES REGISTER 

OLD_BP DU ? 'SAVED BP REGISTER 

OLD.DS DU ? ? SAVED DS REGISTER 

RETURN DD ? 5 RETURN ADDRESS (OFFSET ON TOP, THEN SEGMENT) 

DADDR DD ? ' ADDRESS OF DATA STRUCTURE (OFFSET, THEN SEGMENT) 

ZV DD ? ? ADDRESS OF ZV (OFFSET, THEN SEGMENT) 

XV DD ? S ADDRESS OF XV (OFFSET, THEN SEGMENT) 

PARAMS ENDS 

Sm*******«W*m****m**M«***««******************************* 

DATA SEGMENT PUBLIC 'DATA' 

STOSAVE DU ? 

DATA ENDS 

;***m*##Mm##***#m*m**m****s***M*****M*M#*****m##****** 

CODE SEGMENT PUBLIC 'CODE' 
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08/30/84 PAGE 3 


LCC OBJ 


LINE 

SOURCE 



0000 


106 

FARCOBE PROC 

FAR 




107 






103 






10? 




0000 


110 

RFUN3L: 



0000 IE 


111 

PUSH 

BS 

5SAVE REGISTERS 

0001 55 


112 

PUSH 

BP 


0002 06 


113 

PUSH 

ES 


0003 SBEC 


114 

MOV 

BP, SP 


0005 BS 

R 

115 

MOV 

AX, BATA 


0003 3EC0 


116 

MOV 

ES, AX 


OOQA C55E12 


117 

LBS 

BX, [BP]. XV 

!GET XV ADDRESS AND SEGMENT 

GOOD 9BB907 


118 

FLB 

DWORD PTR EBX] 

.'PUT XV ON TOP OF 8037 STACK 

0010 C55E0E 


119 

LBS 

BX, [BP].ZV 

5GET ZV ADDRESS AND SEGMENT 

0013 9BD907 


120 

FLB 

DWORD PTR [BX] 

.'PUT ZV ON TOP OF 8037 STACK 

0016 C55E0A 


121 

LBS 

BX, [BPLDADBR 

;GET START ADDRESS OF DATA STRUCTURE 

001? 8BB3 


122 

MOV 

BX, BX 

.SAVE BX 

001B 8B2F 


123 

MOV 

BP, EBX3.XPTR 

!B? = XPTR 

001 B 8B4F04 


124 

MOV 

CX, [BX3.MXPTS 

.'CX = NXPTS 

0020 3B7F02 


125 

MOV 

BI, [BX3.ZPTR 

■DI = ZPTR 

0023 D1E7 


126 

SHL 

DI, 1 


0025 D1E7 


127 

SHL 

BI, 1 

5DI = BYTE INDEX = ELEMENT INDEX » 4 

0027 SB7706 


128 

MOV 

SI, CBX3.NZPTS 

iSI = NZPTS 

002A B1E6 


12? 

SHL 

SI, 1 


002C D1E6 


130 

SHL 

SI, 1 

;si = NO. OF BYTES = NZPTS * 4 



131 




002E 


132 

GET.ZJHDEX: 



002E 83C303 


133 

ABB 

BX, 8 

!BX = START OF Z ARRAY 

0031 9BB811 


134 +2 

FCOM 

DWORD PTR CBX + BI] 

; COMPARE ARRAY ELEMENT WITH ZV 

0034 9B26BB3E0000 

R 

135 +1 

FSTSW 

ES: STHSAVE 


003A 9B 


136 +1 

FHAIT 



003B 26A10000 

R 

137 +1 

MOV 

AX, ES: STHSAVE 


003F 250041 


138 +1 

m 

AX, 4100H 

; MASK-IN COMPARISON BITS 

00*2 3D0001 


13? +2 

CMP 

AX, 0100H 


0045 7434 


HO +2 

JE 

LEFTJJF-Z-START 

5 EXIT IF ZV < ARRAY ELE 




HENT 





141 +2 






142 



IF Z > DATAIGRIGINAL INDEX) 

0047 


143 

RIGHT J3F-Z_START: 

*** THEN BEGIN *»* 

0047 83C704 


144 

ABB 

BI, 4 

INDEX = INDEX + 4 

004A 3BFE 


145 

RZSis CMP 

BI, SI 

WHILE (INDEX <= NO. Cf BYTES) 

0040 741B 


146 

JE 

RZS2 




147 



*** DO BEGIN #** 

004E 9BD811 


148 +2 

FCOM 

DWORD PTR [BX + BI3 

COMPARE ARRAY ELEMENT WITH ZV 

0051 9B26DD3EOOOO 

R 

14? +1 

FSTSH 

ES: STHSAVE 


0057 9B 


150 +1 

FUAIT 



0053 26A10000 

R 

151 +1 

MOV 

AX, ES: STHSAVE 


0050 250041 


152 +1 

AND 

AX, 4100H 

: MASK-IN COMPARISON BITS 

005F 3D0001 


153 +2 

CMP 

AX, 0100H 


0062 7405 


154 +2 

JE 

RZS2 

! EXIT IF ZV < ARRAY ELEMENT 



155 +2 




0064 33C704 


156 

ABB 

BI, 4 

INDEX = INDEX + 4 

0067 EB El 


157 

JMP 

RZS1 

: *** END (WHILE) *#* 

0069 8BC6 


153 

RZS2: MOV 

AX, SI 


006B 00 


15? 

SUB 

AX, 4 
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08/30/84 PAGE 4 


LOG OBJ 


LINE 

SOURCE 




006E 83EF04 


160 


SUB 

DI, 4 


0071 3BF9 


161 


CMP 

BIi AX 

IF INDEX O NO. OF BYTES 

0073 7542 


162 


JNE 

FOUNB-LOi’LZ. INDEX 

THEN FOUND J_OW_ INDEX 

0075 9BBBD3 


163 


FSTP 

ST(O) 

ELSE QUT-Cf .RANGE 

0078 EB4B90 


164 


JMP 

GET-X-INDEX 




165 




$** END (IF) *** 



166 




IF Z < DATA( ORIGINAL INDEX) 

007B 


167 

LEFT.OF.Z.START: 

*** THEN BEGIN *** 

007B 83EF04 


163 


SUB 

BI, 4 

INDEX = INDEX - 4 

007E 33FFOO 


169 

LZSis 

CMP 

BI, 0 

WHILE (INDEX > 0) 

0031 7C20 


170 


JL 

LZS2 




171 




«* DO BEGIN 

0033 9BD311 


172 +2 


FCOM 

DWORD PTR tBX + DI3 

COMPARE ARRAY ELEMENT WITH ZV 

0036 9B26BD3E0000 

R 

173 +1 


FSTSW 

ESsSTHSAVE 


008C 9B 


174 +1 


FWAIT 



0030 26ft 10000 

R 

175 +1 


MOV 

AX, ESsSTHSAVE 


0091 250041 


176 +1 


and 

AX, 4100H 

MASK-IN COMPARISON BITS 

0094 3DOOOO 


177 +2 


CMP 

AX, OOOOH 

EXIT IF ZV > ARRAY ELEMENT 

0097 740ft 


178 +2 


JE 

LZS2 


0099 3D0040 


17? +2 


CMP 

AX, 4000H 


009C 7405 


180 +2 


JE 

LZS2 

EXIT IF ZV = ARRAY ELEMENT 

009E 83EF04 


131 


SUB 

BI, 4 

INDEX = INDEX - 4 

00 A 1 EBDB 


182 


Jjsp 

LZS1 

*** END {WHILE) *** 

00ft3 83FF00 


133 

LZS2: 

CMP 

BI, 0 

IF INDEX O 0 

00A6 7B0F 


134 


JGE 

FGUNBJ_OW_Z_ I NDEX 

THEN FOUND-LOW— INDEX 

00.08 83C704 


185 


ABB 

BI, 4 

ELSE OUT-OF-RANGE 

OOAB 9BBDD3 


186 


FSTP 

ST(0> 


OOAE 9BD9EE 


187 


FLBZ 


REPLACE TOP WITH O' 

OOBi 9PD9C9 


183 


FXCH 

SKI) 


00B4 EB1190 


189 


JMP 

GET-X-INBEX 




190 




#« END(IF) *** 



191 





OOB7 


192 

FOUND. 

LOW-Z-INDEX: 


00B7 9BB901 


193 


FLB 

DWORD PTR CBX+BI] 

' STORE ZL ON STACK TOP 

OOBfl 9BBCE9 


194 


FSUB 

ST ( 1 ) , ST 

■(ZV-ZL) 2ND FROM TOP 

OOBB 9BD36904 


195 


FSUBR 

DWORD PTR [BX+DI+41 

; (ZH-ZL) ON STACK TOP 

00C1 9BDEF9 


196 


FBIV 


;{ ZV-ZL )/(ZH-ZL) LEFT ON STACK 

00C4 9BB9C9 


197 


FXCH 

ST ( 1 ) 

?PUT XV ON TOP OF 8037 STACK 



193 





00C7 


19? 

GET— X_ 

INDEX: 



0007 03DE 


200 


ADD 

BX, SI 

;BX = START OF XARRAY 

00C9 8BC7 


201 


MOV 

AX, BI 

;AX = NO. OF Z BYTES 

OOCB F6E1 


202 


MUL 

CL 

;AX = XC ROW OFFSET=N. OF ZBYTES*NXPTS 

OOCB 03B8 


203 


ABB 

BX, AX 

;BX = XC ROW PTR 

OOCF 87FD 


204 


XCHG 

DI, BP 

;BI = XPTR : BP = LOW.ZJNDEX 

GOBI D1E7 


205 


SHL 

DI, 1 


00B3 B1E7 


206 


SHL 

DI, 1 

?DI = BYTE INDEX = XPTR * 4 

00B5 87F1 


207 


XCHG 

SI, CX 

5SI = NXPTS : CX = NO. OF Z BYTES 

00D7 D1E6 


203 


SHL 

SI, 1 

" 

00D9 B1E6 


209 


SHL 

SI, 1 

!SI = NO. OF X BYTES = NXPTS * 4 

OOBB 9BB901 


210 


FLD 

DWORD PTR [BX + DI] 

! STORE XC m TOP of STACK 

CODE 8BC5 


211 


MOV 

AX, BP 


OOEO 050400 


212 


ADD 

AX, 4 


00E3 3BC1 


213 


CMP 

AX, CX 


OOE5 7410 


214 


JE 

GX1 

;exit if z tog high 
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LOG* OBJ 


LINE 

SOURCE 




00E7 03DF 


215 


ABB 

BXi DI 


OOE? 9BD900 


216 


FLB 

DWORD PTR CBX + SIJ 

STORE XB ON TOP OF STACK 

OOEC 2BDF 


217 


SUB 

BX, DI 


OOEE 9BBSE1 


218 


FSUB 

ST, ST(1) 

XB-XC ON ST. TOP 

00F1 9BD3CB 


219 


FHUL 

ST, ST(3) 

<(ZV-ZL)/(ZH-Zl»(XB-XC) ON ST. TOP 

OOF4 9BBEC1 


220 


FABB 


XF ON ST. TOP 

00F7 9BD3D1 


221 +2 

GX1 : 

FCQH 

ST(1) 

COMPARE ARRAY ELEMENT WITH XV 

OOFA 9B26DD3E0000 

R 

222 +1 


FSTSH 

ES:STWSAVE 


0100 9B 


223 +1 


FWAIT 



0101 26A10000 

R 

224 +1 


MOV 

AX, ES:STHSAVE 


0105 250041 


225 +1 


AND 

AX, 4100H 

MASK-IN COMPARISON BITS 

0103 3DOOOO 


226 +2 


CMP 

AX, OOOOH 

EXIT IF XV > ARRAY ELEMENT 

01 OB 74-5C 


227 +2 


JE 

LEFT_OF.X.START 




228 +2 







229 




IF X > DATA (ORIGINAL INDEX) 

010B 


230 

RIGHT.Cf.X.STARTs 

*** THEN BEGIN *** 

010D 83C704 


231 


ABB 

DI. 4 

INDEX = INDEX + 4 

0110 3BFE 


232 

RXSl: 

CMP 

DI, SI 

WHILE (INDEX <= NO. OF BYTES) 

0112 743D 


233 


JE 

RXS2 




234 




*#* DO BEGIN «* 

0114 9BD901 


235 


FLD 

DWORD PTR CBX + DI] 

STORE XD ON TCP OF STACK 

0117 3BC5 


236 


MOV 

AX, BP 


0119 050400 


237 


ABB 

AX, 4 


OilC 3BC1 


233 


CMP 

AX, CX 


01 IE 7410 


239 


JE 

RXS3 

EXIT IF Z TOO HIGH 

0120 03DF 


240 


ABB 

BX, BI 


0122 9BD900 


241 


FLD 

DWORD PTR CBX + SI 3 

STORE XE ON TOP OF STACK 

0125 2BBF 


242 


SUB 

BX, DI 


0127 9BDSE1 


243 


FSUB 

ST, ST ( 1 ) 

XE-XD ON ST. TCP 

012A 9BD3CC 


244 


FHUL 

ST, ST(4> 

( ( Z V-ZL ) / ( ZH-ZL ) ) ( XD— XE ) ON ST. TOP 

012D 9BDEC1 


245 


FABB 


XG ON ST. TOP 

0130 9BD3B2 


246 +2 

RXS3: 

FCOM 

STC2) 

COMPARE ARRAY ELEMENT WITH XV 

0133 9B26BB3E0000 

R 

247 +1 


PSTSK 

ESiSTWSAVE 


0139 9B 


243 +1 


FWAIT 



013A 26A10000 

R 

249 +1 


MOV 

AX, ESsSTHSAVE 


01 3E 250041 


250 +1 


AND 

AX, 4100H 

: MASK-IN COMPARISON BITS 

0141 3DOOOO 


251 +2 


CMP 

AX, OOOOH 

EXIT IF XV > ARRAY ELEMENT 

0144 740B 


252 +2 


JE 

RXS2 




253 +2 





0146 9BB9C? 


254 


FXCH 

ST ( 1 ) 

XF = XG 

0149 9BDDD3 


255 


FSTP 

ST(O) 

' DISCARD OLD XF 

01 4C 83C704 


256 


ADD 

DI, 4 

INDEX = INDEX + 4 

014F EBBF 


257 


JMP 

RXSl 

*** END (WHILE) «* 

0151 9BD9C9 


253 

RXS2: 

FXCH 

ST(l) 


0154 3BC6 


259 


MOV 

AX, SI 


0156 200400 


260 


SUB 

AX, 4 ; 


0159 33EF04 


261 


SUB 

DI, 4 ; 


015C 3BF8 


262 


CMP 

DI, AX 5 

IF INDEX O NO. OF BYTES 

015E 7567 


263 


JNE 

FOUND J.Q&LX- 1 NBEX ; 

THEN FOUND-LON-INBEX 

0160 9BDDDS 


264 


FSTP 

ST (0) 5 

ELSE OUT.Cf .RANGE 

0163 9BBBD3 


265 


FSTP 

ST(O) 

DISCARD XG; AND XV 

0166 EB7490 


266 


JMP 

GETJJNSKER 5 




267 



? 

*** END (IF) «* 



268 



J 

IF X < DATA (ORIGINAL INDEX) 

0169 


269 

LEFOF_X_START: ; 

*** THEN BEGIN *** 
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8086/87/88/186 MACRO ASSEMBLER RFUN3L 


08/30/84 FAOE 


LOC OBJ 

LIME 

SOURCE 




016? 83EF04 

270 


SUB 

DI, 4 5 

INDEX = INDEX - 4 

016C 83FF00 

271 

LXS1 

CMP 

DI, 0 ? 

WHILE (INDEX > 0) 

016F 7C42 

272 


JL 

LXS2 ; 



273 



' 

*** DO BEGIN «* 

0171 9BD901 

274 


FLD 

DWORD PTR [BX + DI) 5 

STORE XD ON TOP OF STACK 

0174 8PC5 

275 


MOV 

AX, BP 


0176 050400 

276 


ADD 

AX, 4 


017? 3BC1 

277 


CMP 

AXi CX 


017B 7410 

278 


JE 

LXS3 

EXIT IF Z TOO HIGH 

017D 03DF 

27? 


ADD 

BX, DI 


017F 9BB900 

280 


FLD 

DWORD PTR CBX + SI] 

STORE XE ON TOP OF STACK 

0182 2BDF 

281 


SUB 

BX, DI 


0184 9BBSE1 

282 


FSUB 

ST, ST(1) 

XB-XC ON ST. TOP 

0187 9BDSCC 

233 


FMUL 

ST, ST(4) 

( ( ZV-ZL ) / ( ZH-ZL > ) ( XB-XC > ON ST. TOP 

018A 9BDEC1 

284 


FADD 

' 

XG ON ST. TOP 

018D 9BD8D2 

285 +2 

LXS3: 

FCOM 

ST ( 2 ) 

COMPARE ARRAY ELEMENT WITH XV 

0190 9B26DD3E0000 R 

286 +1 


FSTSM 

ES:STHSAVE 


0196 ?B 

287 +1 


FWAIT 



0197 26A10000 R 

283 +1 


MOV 

AX, ESiSTUSAVE 


019B 250041 

23? +1 


AND 

AX, 4100H 

MASK-IN COMPARISON BITS 

01 ?E 3D0001 

290 +2 


CMP 

AX, 0100H 


01A1 7410 

2?1 +2 


JE 

LXS2 

EXIT IF XV < ARRAY ELEMENT 

01A3 3B0040 

292 +2 


CMP 

AX, 4000H 


01A6 74 OB 

293 +2 


JE 

LXS2 

EXIT IF XV = ARRAY ELEMENT 

01 A 8 9BD9C9 

294 


FXCH 

ST { 1 ) 

XF=XG 

01AB 9BDDD8 

295 


FSTP 

ST(O) 

DISCARD OLD XF 

01AE 83EF04 

296 


SUB 

DI, 4 

INDEX = INDEX - 4 

01B1 EBB? 

297 


JMP 

LXS1 

*** END (WHILE) *** 

01 B3 83FF00 

293 

LXS2: 

CMP 

DI, 0 

IF INDEX O 0 

01B6 7D0F 

299 


JOE 

FOUND-LCBU-INBEX 

THEN FOUND_LOW-INDEX 

01B8 830704 

300 


ADD 

DI, 4 

ELSE OUT-OF-RANGE 

01BB 9BBBB3 

301 


FSTP 

ST(O) 


01 BE 9BBBB3 

302 


FSTP 

ST (0> 

DISCARD XC 

01C1 9BD9EE 

303 


FLDZ 


REPLACE TOP WITH 0 

01C4 EB1690 

304 


JMP 

GET-ANSWER 



305 




*** END C IF) *** 


306 





01C7 

307 

F0UNB_L01J_X_ INDEX : 

?XF, XG, XV, Z SLOPE ON ST. 

01C7 9BDCEA 

303 


FSUB 

ST(2), ST 

*XF, XG, (XV-XF), ZSLOPE 

01 CA 9BDEE? 

30? 


FSUB 


;(XG-XF), (XV-XF), ZSLOPE 

01 CD 9BDEF? 

310 


FDIV 


;( (XV-XF) /(XG-XFl), ZSLOPE 

01 DO 8BC5 

311 


MOV 

AX, BP 

;IF Z NOT TOO HIGH 

01D2 050400 

312 


ADD 

AX, 4 


01B5 3BC1 

313 


CMP 

AX, CX 

• 

01B7 7403 

314 


JE 

GET-ANSWER 


01D? 9BD9C9 

315 


FXCH 

ST(1) 

! THEN (ZV-ZL) /(ZH-ZL) ON ST. TOP 


316 





01DC 

317 

GET-ANSWER; 



01 DC 8 BC 6 

318 


MOV 

AX, SI 

’AX = NO. OF X BYTES 

01 DE D1E8 

31? 


SHR 

AX, 1 


01 EO D1E8 

320 


SHR 

AX, 1 

? AX = NXPTS = NO. OF X BYTES / 4 


321 




iCX = NO OF Z BYTES 

01E2 F6E1 

322 


MIL 

CL 

'AX = YC ROH OFFSET 


•323 




,’THE AVOVE IN3TR. LIMITS LENGTH OF ZARRAY 


324 




? TO 255 ELEMENTS 
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LOC OBJ 

LIf€ 

SOURCE 



01E4 0303 

325 

ADD 

BX* AX 

iBX = YC RON PTR = NXPTS * ZPTR 

01E6 03DF 

326 

ADD 

BX, BI 

!BX = YC POINTER 

01 ES 9BD907 

327 

FLD 

DWORD PTR EBX] 

! STORE YC ON ST TOP 

01 EB SBG5 

323 

MOV 

AX, BP 


01 ED 050400 

329 

ADD 

AX, 4 


01F0 3BC1 

330 

CHP 

AX, CX 


01F2 7510 

331 

JNE 

Z-NOOIGH 



332 




01F4 

333 

Z.TOQ.HIGH: 


? IF Z OUT OF RANGE 

01F4 8BC7 

334 

MOV 

AX, BI 

•*** DO BEGIN *** 

01F6 050400 

335 

ADD 

AX, 4 

! IF X OUT OF RANGE 

01F9 3BC6 

336 

cmp 

AX, SI 


01 FB 7444 

337 

JE 

FOUND. YV 

! THEN EXIT 

01FD 

333 

X.VALIB: 


! CZ TOO HIGH BUT X IS IN RANGE) 

01FD 9BD94704 

339 

FLD 

DWORD PTR [BX + 4] 

! STORE XD 

0201 EB3290 

340 

JMP 

GET.YV 

>#** END CIF) **# 


341 




0204 

342 

Z.NOT.HIGH: 


; IF Z VALID 

0204 SBC7 

343 

MOV 

AX, DI 

;*** DO BEGIN *** 

0206 050400 

344 

ADD 

AX, 4 

! IF X OUT OF RANGE 

0209 3BC6 

345 

CMP 

AX, SI 

, 

020B 7423 

346 

JE 

Z.VALID 

i THEN EXIT 

02 OD 

347 

BOTH-VALID: 


! ELSE (BOTH XV AND ZV ARE IN RANGE) 

020D 9BD9C9 

343 

FXCH 

ST(l) 

!(ZV-ZL)/(ZH-ZL) ON TOP, THEN YC 

0210 9BD900 

349 

FLD 

DWORD PTR [BX + SI] 

; STORE YB ON ST. TOP 

0213 9BD8E2 

350 

FStlB 

ST, ST(2) 

! (YB-YC) ON ST. TOP 

0216 9BD3C9 

351 

FMUL 

ST, ST ( 1 > 

5 ( ( ZV— ZL ) / ( Zh'-ZL / ) ( YB-YC ) ON ST TOP 

0219 9BDEC2 

352 

FADDP 

ST(2), ST 

iYF 2ND FROM TOP OF STACK 

021C 9B094704 

353 

FLD 

DWORD PTR [BX + 4] 

? STORE YD ON STACK TOP 

0220 9BB94004 

354 

FLD 

DWORD PTR [BX + SI + 43 

; STORE YE ON STACK TOP 

0224 9BD3E1 

355 

FSUB 

ST, ST ( 1 ) 

* YE-YD ON ST TOP 

0227 9BDECA 

•356 

FMULP 

ST (2) , ST 

i(!ZV-ZL)/(ZH-ZL))(YE-YD) 2ND FROM TOP 

022A 9BBEC1 

357 

FADD 


;YG ON TOP OF 8037 

022D EB0690 

353 

359 

JHP 

GET.YV 

;*** end CIF) *** 

0230 

360 

Z.VALID: 


! CX OUT OF RANGE BUT Z VALID) 

0230 03BE 

361 

ADD 

BX, SI 


0232 9BD907 

362 

FLD 

DWORD PTR [BX3 



363 




0235 

364 

GET.YV: 


!YG TOP, THEN YF, THEN SLOPE 

0235 9BDSE1 

365 

FSUB 

ST, ST ( 1 ) 

JYG-YF ON ST. TCP 

0238 9BDECA 

366 

FMULP 

ST(2), ST 

; ( ( X V-XF ) / ( XG-XF ) ) ( YG-YF ) 2ND FROM TOP 

023B 9BBEC1 

367 

FADD 


!Y V ON ST. TOP 

023E EB0190 

368 

JMP 

FOUfflLYV 



369 





370 




0241 

371 

FOUND. W: 



0241 

372 

SAVEJOLPOINTERS: 


0241 8BDA 

373 

MOO 

BX, DX 

i RESTORE BX TO START OF DATA STRUC. 

0243 D1EF 

374 

SHR 

DI, 1 

0245 D1EF 

375 

SHR 

DI, 1 


0247 S93F 

376 

MOV 

[BXLXPTR, DI 

ISAVE XPTR 

0249 DIED 

377 

SHR 

BP, 1 


024B DIED 

373 

SHR 

EP, 1 


024D 896F02 

379 

MOV 

[BXLZPTR, BP 

iSAVE ZPTR 
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LOC OBJ LINE SOURCE 


0250 

0250 07 

0251 5D 

0252 IF 

0253 CAOCOO 


380 

381 EPILOGUE: 

882 POP ES 

333 POP BP ; RESTORE REGISTERS 

384 POP DS 

335 RET 12 

386 

387 FARCOBE ENDP 
333 

38? CODE ENDS 

390 

393 

394 END 


ASSEMBLY COMPLETE, NO ERRORS FOUND 


48 





B 

Figure 5. - Map function interpolation diagram. 


16-Bit segment register 

Regular 16-bit register 
1 1 0 1 1 0 1 1 0 0 0 1 1 1 1 0 






FIELD 


DEFINITION 


XPTR: 

ZPTR: 

NXPTS: 

NZPTS: 

Data Arrays: 


X pointer 
Z pointer 

Number of x data points per z curve 
Number of z curves 

These fields hold actual experimental data points. 
The points mark the start and end of a curve 
segment. 


Figure 7. - Field identifications in breakpoint information block. 


DIMENSION XARRAY (7). YARRAY (7) 

INTEGER XPTR, YPTR 

COMMON /STUFF/XPTR, ZPTR, NXPTS, NZPTS, XARRAY, YARRAY 
DATA XPTR, ZPTR, NXPTS, NZPTS/0,0, 7, 7/ 

DATA XARRAY /0. 0, 1. 0, 2. 0, 3. 0, 4. 0, 5. 0, 6. 0/ 

DATA YARRAY /0.0, 1.0, 4.0, 9.0, 16.0, 25.0, 36.0/ 

XIN- 1.5 

Y- RFUN1L (XIN, XPTR) 

PRINT* Y 
END 

Figure 8. - Fortran program for univariate function call. 


y 



Figure 9. - Two out-of-range examples. 



FIELD NAME FIELD VALUE 


XPTR 

X index 

ZPTR 

Not used 

NXPTS 

Number of x data points 
(4 in this example) 

NZPTS 

Not used 

XARRAY 

XI, X2, X3, X4 

YARRAY 

Yl. Y2, Y3, Y4 


Figure 10. - Breakpoint information block for RFUN1 programs. 


FIELD NAME 

FIELD VALUE 

XPTR 

X index 

ZPTR 

Z index 

NXPTS 

Number of x data points 
( 4 in this example) 

NZPTS 

Number of z data points 
(5 in this example) 

ZARRAY 

Zl, Z2, Z3, Z4, Z5 

XARRAY 

XI, X2, X3, X4 

YARRAY 

Yll. Y12, Y13, Y14 
Y21, Y22, Y23, Y24 
Y31, Y32, Y33, Y34 
Y41, Y42, Y43, Y44 
Y51, Y52, Y53, Y54 

Figure 11. - 
grams. 

Breakpoint information block for RFUN2 pro- 


FIELD NAME 

FIELD VALUE 

XPTR 

X index 

ZPTR 

Z index 

NXPTS 

Number of x data points 
(4 in this example) 

NZPTS 

Number of z data points 
(5 in this example) 

ZARRAY 

Zl, Z2, Z3, Z4, Z5 

XARRAY 

Xll, X12, X13, XI 4 
X21, X22, X23, X24 
X31, X32, X33, X34 
X41, X42, X43, X44 
X51, X52, X53, X54 

YARRAY 

Yll. Y12, Y13, Y14 
Y21, Y22, Y23, Y24 
Y31, Y32, Y33, Y34 
Y41, Y42, Y43, Y44 
Y51, Y52, Y53, Y54 


Figure 12. - Breakpoint information block for RFUN3 
programs. 
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